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Executive  Summary 


In  response  to  a  recommendation  stemming  from  an  earlier  report*  that 
“the  knife  edge  model  be  extended  to  include  multiple  knife  edges,”  this 
study  develops,  as  a  special  case,  a  computer  program  model  of 
microwave  diffraction  over  two  edges/hills  that  are  situated  on  a  flat  earth. 


In  summary,  the  following  has  been  accomplished  in  this  report: 

An  extension  is  made  of  the  Meeks  knife-edge  method  for  computing 
power  density  on  target,  in  order  to  consider  double  diffraction.  Whereas 
Meeks’  method  (which  itself  is  an  extension  of  Deygout’s  1966 
simplification  of  diffraction  computations  for  the  Kirchhoff  knife-edge 
theory)  models  the  broadcast  of  microwaves  over  a  single  knife  edge,  the 
method  described  here  models,  at  the  discretion  of  the  user,  either  isotropi¬ 
cally  broadcast  or  anisotropically  beamed  microwave  radiation  over  two 
knife  edges. 


In  order  to  make  this  extension,  one  constructs  a  set  of  uniquely  described 
ray  paths  from  the  source  antenna  to  the  sink  antenna.  The  model  de¬ 
scribes  double  diffraction  of  microwaves  across  electrically  interactive  ter¬ 
rain  with  ground  reflections  at  low  incidence  angles.  This  method  can  be 
further  extended  for  diffraction  over  more  than  two  edges. 


* Albert  G.  Gluckman,  The  Lobing  Structure  of  ML  • owave  Radiations  due  to  Reflection  and  Diffraction  from 
Terrain  (U),  Harry  Diamond  Laboratories,  HDL-TM-86-10  (September  1986).  (CONFIDENTIAL) 
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1.  Introduction 


Suppose  a  ray  path  of  a  microwave  train  is  diffracted  by  two  knife  edges 
which  are  set  on  a  flat  earth.  How  can  we  determine  by  how  much  or  by 
how  many  decibels  the  power  level  is  diminished  by  the  diffraction?  And 
if  interference  occurs  between  this  ray  and  a  direct  ray,  how  does  this  in¬ 
terference  affect  the  power  level  of  the  train  of  microwaves?  Further,  if 
reflection  occurs  at  one,  two,  or  three  points  that  are  separated  by  the  knife 
edges  (see  fig.  1),  how  much  power  density  of  the  train  of  microwaves  is 
lost  to  eddy  currents  in  the  conductive  soil? 

Apropos  both  the  computational  characteristics  and  the  theoretical  under¬ 
pinnings/foundations  of  propagation  models,  the  following  question  comes 
to  mind.  How  can  the  reduction  in  power  level  for  a  single  knife  edge 
[1-6]*  be  extended  to  apply  to  the  problem  of  determining  the  reduction  in 
power  level  for  two  knife  edges,  for  all  reflection  path  geometries  that  are 
encountered? 

The  purpose  of  this  report  is  to  show  how  this  may  be  accomplished.  The 
FORTRAN  program  to  do  this  (called  DOUBLEDGE)  is  shown  in  appen¬ 
dix  B.  It  can  handle  the  propagation  of  either  broadcast  radiation  at  vhf  or 
of  beamed  microwaves  which  have  much  higher  frequencies. 


Figure  1.  Knife-edge 
diffraction  and  re¬ 
flection  for  example 
of  two  diffractions 
and  three  reflections. 


Target 


I:  1st  point  of  reflection  on  terrain 
IT-  2d  point  of  reflection  on  terrain 
III:  3rd  point  of  reflection  on  terrain 
Ff.  Edge  of  1st  diffraction 
^2  ■'  Edge  of  2d  diffraction 


* References  are  listed  at  the  end  of  the  main  body  of  text  (p  23). 
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The  report  first  discusses  propagation  of  the  direct  ray,  next  the  indirect 
ray,  and  then  the  application  of  the  double-knife-edge  model,  ending  with 
a  discussion  of  the  results. 

Appendix  A  shows  the  diffractive  ray-path  geometries  for  two  knife  edges, 
appendix  B  is  a  listing  of  the  FORTRAN  computer  program  of  the  model, 
and  appendix  C  is  a  program  listing  which  was  used  to  extract  the  data 
subset  for  making  the  plots  shown  in  appendix  D. 

The  plots  in  appendix  D  show  the  change  in  power  density  at  the  target, 
expressed  in  decibel  scale  (relative  to  free  space  level  of  0-dB  decrement), 
versus  target  height.  The  plots  show  a  good  agreement  of  accepted  theory 
with  the  results  from  the  model.  Because  of  its  simplicity,  this  method 
takes  little  computer  time. 

Appendix  E  contains  plots  of  single  diffraction  and  double  diffraction  of 
microwaves.  Some  of  these  plots  are  used  to  show  how  power  density  at 
the  target  after  single  diffraction  will  compare  to  power  density  at  the  tar¬ 
get  after  double  diffraction.  Some  other  plots  in  appendix  E  show  how,  by 
changing  emitter  height,  battlefield  saturation  by  radiation  can  be 
enhanced. 

Appendix  F  shows  a  plot  comparing  a  theoretically  derived  curve  for  two 
diffractions  over  electrically  interactive  terrain  with  a  curve  that  is  derived 
from  actual  propagation  data  collected  in  the  field  from  measurements. 
These  propagation  data  are  taken  with  respect  to  the  transmitter  on  an 
azimuth  to  a  location  called  Forest  Hill. 

Appendix  G  shows  a  plot  comparing  a  theoretically  derived  curve  for  a 
single  diffraction  over  terrain  (using  the  KNIFEDGE  program)  to  a  curve 
prepared  from  field  measurement  data  for  two  diffractions  over  terrain. 
These  propagation  data  are  taken  with  respect  to  the  transmitter  on  an 
azimuth  to  Forest  Hill. 


2.  Direct  Ray  Path 

There  are  three  possible  geometries  to  describe  the  clearance  of  the  direct 
ray  path  over  two  knife  edges.  The  three  geometries  depend  on  whether 
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the  heights  of  the  two  knife  edges,  Fl  and  F2,  are  as  F}  <  F2,  as  F}  >  F2,  or 
as  F j  =  F2-  These  geometries  are  pictured  in  figure  2. 

With  sufficient  distance  above  the  knife  edges,  relative  to  wave  length, 
these  ray  paths  exhibit  no  diffraction.  They  are  called  the  direct  ray  paths, 
or  the  direct  rays. 

Let  us  consider  knife-edge  heights  Fx  and  F2;  the  distances  dv  d2,  and  d3 
downrange;  the  height  of  the  transmitter,  zx;  and  the  height  of  the  target, 
z2,  as  illustrated  in  figure  2.  From  these  three  diagrams  in  figure  2,  a  logi¬ 
cal  choice  can  be  made  as  to  which  path  geometry  is  appropriate.  This  is 
done  by  a  testing  algorithm  which  is  based  on  (1)  the  sign  and  magnitude 
of  the  slope  angles  of  the  paths  over  the  edges  and  on  (2)  whether  F,  is 
greater  than,  equal  to,  or  less  than  F2. 

From  figure  3,  the  slope  angles  are 

ai  =  tan-1  [(Fi  -zi)/di] 


and 


So,  if 


ct2  =tan_1  [(F2  -  zi)/(d\  +  d2)]  . 


Fj  <  F2  and  ctj  <  ctj  , 

the  direct  path  has  F2  as  its  edge  of  closest  approach.  In  this  case,  the  nor¬ 
malized  electric  amplitude  A 12  associated  with  the  Fresnel  diffraction  of  an 
assumed  wave  over  F2  is  calculated.  If 

Fj  <  F2  and  (Xj  2:  , 


or  if, 


F j  >  F2  and  otj  t  , 

the  direct  path  has  Fj  as  its  edge  of  closest  approach.  In  this  case,  the  nor¬ 
malized  electric  amplitude  An  associated  with  the  Fresnel  diffraction  of  an 
assumed  wave  over  Fx  is  calculated. 
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Figure  2.  Clearance  of  direct  ray  path  over  knife  edges  F,  and  F. 


F,  >  F2 
ot,  >  a2 


Figure  3.  Geometry  for  determining  direct  path  over  edge  F,  or  Fv  based  on  slope  angle  and  ray 
path  and  relative  heights  of  Fx  and  Fr 


The  normalized  electric  amplitudes,  A,  of  these  direct  ray  paths  (as  well  as 
the  indirect  paths  discussed  in  the  next  section,  which  experience  reflec¬ 
tion  as  well  as  diffraction)  are  computed  with  the  aid  of  the  Cornu  spiral 
approximations  for  the  Fresnel  integrals. 


3.  Indirect  Ray  Path 

The  indirect  ray  path  is  a  path  along  which  a  microray  experiences  dif¬ 
fraction  from  one  or  more  knife  edges,  as  well  as  reflection  from  the 
grcund/terrain.  A  microray  is  a  ray  at  microwave  frequency.  The  theory 
of  ray  optics  was  applied  by  Deygout,  and  later  Meeks,  to  the  problem  of 
microwave  propagation  from  a  point-source  emitter. 

The  ray  path  emanating  from  the  point-source  microwave  transmitter  is 
hypothesized  to  be  a  continuous  signal.  This  signal  can  be  diffracted  by 
intervening  hills/ridges,  as  well  as  reflected  by  flattened  terrain  profile  fea¬ 
tures.  Interference  patterns  arise  when  these  diffracted  and  reflected  rays 
interfere  with  the  direct  rays. 

This  method  for  determining  diffraction  over  two  geometries  takes  into  ac¬ 
count  the  diminution  of  electric  amplitude  that  is  associated  with  each  ray 
path  geometry. 

From  figures  2  and  4,  it  can  be  seen  that  reflection  of  a  ray  can  occur 

on  dl  or  on  d2  or  on  d3  separately; 
or  on  dx  and  d2,  or  on  dx  and  d3,  or  on  d^  and  d3  conjointly; 
or  on  dx  and  d2  and  d^  in  concert. 

Furthermore,  it  is  evident  that  diffraction  can  occur  at  F }  or  at  F2  or  at  both 
F j  and  F2,  together  with  any  of  the  possibilities  of  reflection  just  men¬ 
tioned.  This  means  that  a  total  of  14  ray-path  geometries  are  possible  for 
an  indirect  ray  path. 

Another  level  of  complication  in  the  theory  is  that  the  target  height  is 
raised  in  discrete  intervals,  so  that  for  each  increase,  the  reflection  points 
alter  or  shift  in  position  on  the  terrain,  and  the  grazing  angles  change. 
Since  the  reflecting  electric  ray  interacts  electrically  with  the  terrain,  eddy 
currents  are  formed  in  the  soil  at  the  point  of  reflection,  and  electric  energy 
is  lost  in  the  soil  in  the  form  of  Joule  heat.  This  means  that  the  resulting 
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reflected  electric  microray  has  a  diminished  amplitude.  The  reflection 
coefficient  depends  upon  permittivity  (which  is  the  dielectric  constant)  of 
the  soil,  wave  length,  conductivity  of  the  soil  (in  siemens),  and  the  grazing 
angle  (in  radian  measure).  Figure  4  shows  how  to  determine  the  re¬ 
spective  grazing  angles  PS/1,  PS/2,  and  PS/3  of  a  ray  over  terrain  regions 
dv  d2,  and  d3,  respectively. 

In  the  computer  program  (with  fig.  3  in  mind),  CTV  is  the  vertical  polar¬ 
ization  reflection  coefficient  and  CTH  is  the  corresponding  horizontal 
polarization  reflection  coefficient.  The  reflection  coefficient  determines 
the  diminution  (or  reduction)  in  the  (normalized)  electric  amplitude  of  the 
ray,  after  reflection  on  the  ground.  This  means  that  after  reflection,  there 
is  a  corresponding  reduction  in  the  eventual  power  density  output  at  target 
from  this  ray  path. 

The  subroutine  FRESNL  (in  app  B)  computes  the  reflection  coefficient. 
The  FORTRAN  parameter  ANG  in  subroutine  FRESNL  takes  the  values 
PS/1,  PS/2,  or  PS/3  of  the  grazing  angles  in  radians.  PS/1,  PS/2,  and  PS/3 
are  the  grazing  angles  over  dv  d2,  and  dv  respectively. 


Figure  4.  Ray-path 
geometry  for  deter¬ 
mining  grazing  an¬ 
gles  over  terrain  re¬ 
gions  (needed  for 
computing  reflection 
coefficients  CTp 
where  i  =  1,  2,  or  3). 


(a)  over  terrain  region  d: 


(b)  over  terrain  region  d2 


-/ 

tan  (z/  +  Fj)/di  =  PSI/ 


-I 

tan  (Z2  +  Fj)/fd:  +  dii  =  PSP 
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For  either  horizontally  or  vertically  polarized  microwaves,  the  calculated 
Fresnel  reflection  coefficient  ( CTH  or  CTV)  and  the  ground  reflectivity 
REFLi  (in  FORTRAN  computer  code)  are  multiplied  together  in  one-to- 
one  correspondence  with  the  appropriate  terrain  interval  dt.  This  means 
that 


CTi  = 


x  REFLi 


where  /  =  1 ,  2,  or  3  , 


and  CT-  is  the  amended  reflection  coefficient. 

In  order  to  distinguish  beamed/anisotropic  propagation  at  microwave  fre¬ 
quency  from  broadcast/isotropic  propagation  at  vhf,  the  computer  program 
relies  on  the  data  input  of  either  “anisotropic”  or  “isotropic.”  If  a  beam  is 
formed,  a  corresponding  high-frequency  input  is  necessary  in  order  to 
satisfy  the  physical  requirements  of  the  generic  theory.  Just  as  in  the 
KNIFEDGE  program  [7]  of  the  diffraction  of  a  microray  (or  its  associated 
microwave)  caused  by  a  single  edge,  so  too,  in  this  DOUBLEDGE  pro¬ 
gram  model  is  the  nonisotnopic  radiator  incorporated  into  the  theory  by 
considering  the  on-axis  ray  of  the  antenna  as  the  direct  ray  and  the  off-axis 
ray  as  the  reflected  ray.  The  magnitude  of  the  reflected  ray  is  less  than 
that  of  the  direct  ray  by  an  E-field  reduction  factor: 

Er  =  GRdB  =  10(G-C)/20  , 


where  G  is  the  gain  of  the  antenna  in  decibels  and  C  is  the  on-axis  antenna 
gain.  In  the  isotropic  case,  gain  due  to  directivity  is  not  considered. 

Corresponding  to  each  of  the  two  knife  edges  Fx  and  F2,  there  is  a  clear¬ 
ance  parameter  of  a  first  Fresnel  zone.  For  edge  F v  this  is 

*  _  /  ^1  ^2 
\  3 TTd2 


and  for  F2,  this  is 


&23 


"Kdj  dy 
di  +  dy 
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For  the  edge  F v  the  argument  u  of  the  Fresnel  integrals 

u  u 

C(u)  =  J*  cos  j  u2  du  and  S(u)  -  J  sin  ^  u2  du 
0  0 

when  multiplied  by  the  factor  2~m  forms  the  ratio  of  the  clearance 
parameter  A  of  the  direct  ray  over  or  through  (if  negative  clearance  below) 
the  knife  edge  Fj  expressed  in  units  of  clearance  of  the  First  Fresnel  zone 
A12.  Thus  m  =  A  21/2/A12. 

On  the  other  hand,  for  the  edge  F2,  u2~1/2  is  the  ratio  of  the  clearance 
parameter  A  of  the  direct  ray  over  or  through  (if  negative  clearance  below) 
the  knife  edge  F2  expressed  in  units  of  clearance  of  the  first  Fresnel  zone. 
Thus  u  =  A  2ml  A^. 

Clearance  parameters  A12,  corresponding  to  edge  Fv  are  calculated  for 
each  of  the  seven  generic  indirect  ray  paths  over  Fv  and  this  calculation  is 
iterated  for  each  change  that  occurs  due  to  the  stepping  of  the  target  height 
as  it  gains  altitude.  Likewise,  clearance  parameters  A23,  corresponding  to 
edge  F2,  are  calculated  for  each  of  the  defined  seven  generic  indirect  ray 
paths  over  F2. 

The  choice  of  the  direct  ray  path  may  be  taken  over  either  knife-edge  F [  or 
F2,  but  this  choice  is  decided  by  the  method  that  is  described  in  section  2. 
If  knife  edge  Fj  is  the  one  of  closest  approach  to  the  ray,  yet  far  away 
enough  not  to  contribute  any  diffractive  effects  (i.e.,  the  bending  of  the  ray 
due  to  the  presence  of  the  knife  edge),  the  normalized  electric  amplitude 
An,  corresponding  to  F x  is  calculated.  If,  on  the  other  hand,  F2  satisfies 
this  condition,  the  normalized  electric  amplitude  A]2  is  calculated.  Note 
that  in  the  notation  An  and  A12,  the  first  subscript,  1,  refers  to  the  direct 
ray;  the  second  subscript  refers  to  either  edge  Fj  or  F2. 
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In  more  detail,  normalized  electric  field  amplitudes  for  each  of  the  direct 
rays  over  F{  or  F2  are  computed  independently;  and  this  is  true  for  all 
computations  involving  the  indirect  rays  as  well.  The  choice  of  either  di¬ 
rect  ray  over  Fl  or  over  F2  is  made  by  the  decision  algorithm  based  on 
slope  and  the  relative  heights  of  edges  F l  and  F2  that  is  described  in  figure 
3  of  section  2.  There  can  only  be  one  direct  ray. 

Likewise,  normalized  electric  field  amplitudes  Atj  corresponding  to  each  of 
the  14  remaining  generic  indirect  ray  paths  are  computed  in  a  stepwise 
fashion  for  each  iteration,  as  the  target  height  increases  step  by  step  to  the 
maximum  designated  height.  The  electric  field  contributions  for  all  paths 
at  a  particular  height  are  summed  for  their  total  field  value.  The  normal¬ 
ized  electric  amplitudes  A  -,  where  i  =  1  through  8  (corresponding  to  path), 
and  j  -  1  or  2  (corresponding  to  edge  F x  or  F2),  is  derived  from  the 
equation 

=  (C  +  0.5)2  +  (S  +  0.5)2  , 

where 

0  ZAij  <  1  . 

All  Ay  values  are  derived  independently  of  each  other.  For  the  case  of  the 
14  paired  indirect-ray-path  geometries,  the  electric  amplitudes  of  the  seven 
pairs  of  ray  paths  (shown  in  app  A)  are  combined  as 

A21A22,  A31A32,  Am  A42 »  A51A52,  i4gi  Af,2  *  An\Ai2, 

for  all  iterations  on  target  height.  These  products  of  normalized  electric 
amplitudes  reduce  the  available  energy  on  target,  because  of  diffractive 
loss  due  to  the  bending  of  rays  into  the  shadow  region  behind  the  hills. 

Because  the  values  of  the  amplitudes  Ay  theoretically  lie  between  0  and  1 
(except  that  in  actual  practical  use  numerical  error  may  cause  this  condi¬ 
tion  to  be  sometimes  violated  to  some  extent),  they  act  to  diminish  the 
power  density  of  the  electric  field  corresponding  to  each  ray  path  associ¬ 
ated  with  this  diffracted  foward  scattering.  The  combined  effect  of  the 
pairwise  products  of  these  amplitudes  (which  may  be  called  product 


amplitudes),  as  shown  just  above,  are  a  measure  of  the  reduction  in  power 
density  at  the  target  due  to  the  diffraction  induced  by  both  knife  edges  Fx 
an  dF2. 

The  expression  CEi  is  a  complex  number  expressing  the  magnitude 
(amplitude)  and  the  phase  lag  angle  of  each  contribution  of  the  electric 
field  corresponding  to  each  iteration  of  the  target  height,  for  each  generic 
path.  The  total  magnitude  for  each  field  contribution  is  expressed  as 
AnAi2  multiplied  by  the  appropriate  combination  of  reflection  coefficients, 
CT-,  and  by  the  gain  reduction  factor,  ER,  to  account  for  the  beam  con¬ 
centration  of  microwave  radiation. 

The  phase  lag  angle  B,  where 

n  *  -1  ( S  +0.5  V  0 

fi=tan  |^CTo3 J  +  T  •  *="lor3  - 

is  due  to  the  difference  between  the  phase  angle  of  the  reflected  ray  and  di¬ 
rect  ray,  at  target,  and  this  is  a  consequence  of  the  differences  in  the  path 
lengths  of  the  direct  and  indirect  (reflected)  rays.  Therefore,  the  relative 
power  density  at  target  in  decibel  units,  relative  to  propagation  in  free 
space  is 

8 

20  log  C  E\  +  S  C  Ei  , 

i=2 

where  CEl  is  the  complex  form  for  the  electric  contribution  Au  or 
A12  e B'f~i  for  the  direct  ray;  and  C£,,  i  =  2  to  8,  is  the  complex  form  of  the 
electric  contribution  for  each  of  the  paired  generic  indirect  ray  paths,  ex¬ 
pressed  as 


CE22  =A2tA22  e8^ e'2nARil,xe-2nARnlxri  Er 
CE2 2  =  A31 A32  e8^ e~2^llXe-2n^32lXri  Er 
CEa2  =  AiiA^ze8^ e-2****i*e-***»t\rir2  Er 
CEsz  =A5lA52eB'lz'e-2K*R>',xe-2nAR>2lXri  Er 
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CE6 2  =  A6iA62eB^ e-^'^e-^^TzTs  ER 

ce12  =  a~!\  a12  e* ^  e-^ixe~2^llixrx  r3  ER 

CES2  =AixAi2eB'rxe-2l**",'K-e-2l*RnlxTxr2Ti  ER 

where  the  amended  reflection  coefficients  CTi  =  T-  each  correspond  to 
their  associated  terrain  region  dt. 

The  normalized  electric  amplitudes  A(1  and  Ai2  ( i  =  2  through  8)  refer  to 
diffraction  at  edges  Fx  and  F2,  respectively,  for  all  pairs  of  ray  path 
geometries  that  are  described  in  figure  A-l  of  appendix  A. 


4.  Application  of  the  DOUBLEDGE  Computer  Program 

Figure  5  shows  how  terrain  with  two  prominent  features/hills  can  be  pic¬ 
tured  for  the  particular  example  where  the  transmitter,  T,  is  set  up  on  a 
high  hill  overlooking  the  predominantly  downhill  terrain.  Two  large  hills 
lie  between  the  transmitter,  T,  and  the  target,  5,  located  downrange.  The 
intervening  hilltops  are  located  at  positions  (X4,  f4)  and  (X6,  f6).  The  line 
XjX2  is  parallel  to  sea  level  in  the  flat  earth  model.  The  terrain  distance 
and  height  parameters  are  referenced  to  the  lowest  position  above  sea 
level,  X2,  in  the  same  manner  as  this  was  done  by  Gluckman  [7]  in  his 
KNIFEDGE  computer  model  of  a  microwave  beam  that  is  diffracting  over 
a  single  hill/ridge.  The  microwave  targeting  position  lies  above  position 
X2  at  location  ground  zero. 

With  the  use  of  the  appropriate  data  derived  from  actual  measured  field 
transmissions  and  receptions,  this  doubly  diffractive  model  with 
attendant/accompanying  reflections  can  be  tested  for  its  accuracy  in  repre¬ 
senting  actual  physical  conditions  of  terrain,  polarity  of  the  wave,  and  fre¬ 
quency,  with  respect  to  the  measured  signal  strength. 
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Figures.  Double' 
knife-edge  model 
over  terrain. 


0C5.  Ys) 


T  at  zj 


5.  Discussion  of  Results 

By  the  expedient  of  raising  or  lowering  the  height  of  the  high-power 
microwave  transmitter,  one  can  alter  the  elevations  of  the  power  density 
lobes.  This  would  result  in  increasing  the  coverage  of  electric  energy  on 
targets  over  or  about  the  position  of  ground  zero.  This  procedure  would 
change  the  amount  of  diffracted  energy  in  the  umbra/shadow  which  is  be¬ 
low  the  line  of  sight,  thereby  increasing  the  energy  level  of  radiation  that 
can  be  placed  on  a  target  behind  a  hill.  The  method  developed  in  this  re¬ 
port  to  determine  power  loss  resulting  from  two  diffractions,  for  ray  paths 
reflecting  from  electrically  interactive  terrain,  can  be  extended  to  the  case 
of  diffraction  occurring  over  n  edges. 

5.1  Comparison  of  Power  Loss  of  Double  Versus  Single 
Diffraction 

In  comparing  the  plots  of  appendix  E  in 
figure  E- 1 ,  where  D 1  =  2500  m  and  D2  =  5000  m, 
figure  E-2,  where  D  \  -  3750  m  and  D2  =  3750  m,  and 
figure  E-3,  where  D 1  =  5000  m  and  D2  =  2500  m, 
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all  other  electric  and  geometric  terrain  parameters  remaining  the  same,  one 
makes  the  following  observations. 

In  all  three  plots,  power  loss  due  to  double  diffraction  is  greater  than  for 
single  diffraction,  and  the  double-diffraction  power  loss  curves  asymptoti¬ 
cally  follow  the  single  diffraction  curves  for  dB  down  (loss).  Above  line 
of  sight  at  target,  the  power  loss  curves  for  double-  and  single-edge  dif¬ 
fraction  intertwine  with  each  other  in  a  dampening  oscillation  through/at 
the  0-dB  region.  It  is  noticed  that  when  the  single  knife  edge  is  closer  to 
the  emitter,  the  power  loss  curves  for  double-  and  single-edge  diffraction 
show  a  closer  asymptotic  approach,  as  shown  in  figure  E-l.  As  the  single 
knife  edge  arrives  at  the  midpoint  of  the  target  range,  the  two  curves  fall 
further  apart,  nevertheless  maintaining  their  parallel  character,  as  shown  in 
figure  E-2.  This  characteristic  is  further  accentuated  as  the  single  knife 
edge  approaches  the  target,  as  shown  in  figure  E-3.  In  all  three  cases, 
however,  the  power  loss  curve  at  target,  of  diffraction  due  to  two  knife 
edges,  exhibits  spikes  of  high  loss,  but  these  spikes  are  very  narrow. 
However,  these  spikes  appear  as  natural  extensions  of  the  oscillations  of 
the  single-knife-edge  power  loss  curves.  The  tails  of  these  isolated  spikes 
may  show  a  power  loss  at  their  corresponding  height  above  the  terrain,  as 
much  as  30  to  35  dB  less  than  power  loss  due  to  single-edge  diffraction. 

5.2  Comparison  of  Power  Loss  of  3-m  Versus  6-m  Antenna 
Height 

Diffraction  due  to  two  knife  edges.  Plot  E-4  (see  app  E)  shows  that  power 
loss  is  greater  for  an  emitting  antenna  placed  at  a  height  of  3  m  than  at 
6  m.  Note  that  in  both  cases,  all  the  electric  and  geometric  terrain 
parameters  remain  the  same. 

Diffraction  due  to  a  single  knife  edge.  In  comparing  the  plots  of  appendix 
E  in 

figure  E-5,  where  D 1  =  2500  m  and  D2  =  5000  m, 
figure  E-6,  where  Dl  =  3750  m,  and  D2  =  3750  m, 
figure  E-7,  where  D 1  =  5000  m  and  D2  =  2500  m, 

all  other  electric  and  geometric  terrain  parameters  remaining  the  same,  one 
makes  the  following  observation. 
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Where  the  knife  edge  is  at  first  closer  to  the  emitter,  and  moves  to  a  loca¬ 
tion  midway  between  the  emitter  and  the  target,  as  shown  in  figures  E-5 
and  E-6,  the  power  loss  is  greater  for  an  antenna  set  at  a  3-m  height  than 
for  an  antenna  situated  higher,  at  a  6-m  height  above  the  terrain.  Both 
power  loss  curves,  however,  show  parallelism,  as  if  they  were  shifted  apart 
by  a  horizontal  translation  along  the  abscissa.  At  the  0-dB  vertical  region 
of  the  plot  (at  free  space  propagation  on  the  decibel  scale  considered  as  a 
reference-value/origin,  against  which  to  measure  propagation  over  terrain), 
the  two  power  loss  curves  intertwine  in  a  dampening  oscillation  as  altitude 
of  the  target  increases  (i.e.,  as  altitude  increases,  the  dampening  of  the 
power  loss  oscillation  increases). 

Where  the  knife  edge  is  set  farther  away  from  the  emitter,  and  therefore 
closer  to  the  target,  as  shown  in  figure  E-7,  the  difference  in  diffractive 
power  loss  due  to  changing  the  emitter-antenna  height  from  3  to  6  m  is 
lost. 

5.3  Verification  of  the  DOUBLEDGE  Diffraction  Theory  with 
Collected  Field  Data 

In  November  1979  and  February  1980,  M.  L.  Meeks  conducted  an  on-site 
measurement  survey  of  electric  field  output  (or  equivalently,  of  power  out¬ 
put)  for  the  Lincoln  Laboratory  of  the  Massachussetts  Institute  of  Technol¬ 
ogy  (MIT).  His  transmitter  was  a  vhf  omnidirectional  radio  station,  used 
as  an  aircraft  navigation  aid,  situated  on  a  hilltop  about  80  km  west  of 
Boston,  MA.  This  station  transmitter  propagated  a  signal  at  1 10.6  MHz, 
with  an  omnidirectional  pattern  which  was  symmetrical  in  azimuth.  The 
polarization  was  horizontal. 

The  survey  measurements  were  made  onboard  a  vertically  descending 
helicopter.  The  receiving  antenna  was  a  horizontal  dipole  located  under 
the  fuselage.  The  effect  of  antenna  gain  due  to  elevation  was  ignored.  In 
the  horizontal  plane,  however,  gain  was  isotropic,  that  is,  the  same  in  all 
directions.  Additional  details  of  these  survey  measurements  are  available 
elsewhere  [4,5]. 

Lincoln  Laboratory  supplied  Harry  Diamond  Laboratories  (HDL)  with 
raw  digital  measurement  data  from  the  azimuthal  propagation  path,  along 
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which  (for  our  study)  the  signal  was  measured  above  the  terrain  at  the  tar¬ 
get.  This  path  was  from  the  transmitter  to  a  location  called  Forest  Hill. 
Two  measurements  were  made  independently  for  this  path  across  the  ter¬ 
rain.  This  particular  path  was  chosen  because  it  displayed  two  prominent 
diffractive  hills  across  which  the  signal  was  propagated. 

The  forest  population  consisted  of  a  mix  of  evergreen  and  deciduous  trees 
which  had  lost  their  leaves  for  the  winter  season. 

At  HDL,  the  two  independent  aggregates  of  signal  measurements  were 
averaged  at  5-m  interval  heights  for  this  path  from  the  transmitter  to  the 
Forest  Hill  location. 

The  curve  shown  in  plot  F-l  of  appendix  F,  describing  an  averaged  data 
measurement  set,  was  compared  with  the  theoretically  derived  curve  for 
this  azimuth.  As  can  be  seen,  the  results  are  very  favorable,  and  the  close¬ 
ness  of  fit  verifies  the  validity  of  the  “doubledge”  theory.  It  is  recom¬ 
mended  that  a  statistical  study  be  made  of  the  characteristics  of  this 
closeness. 

In  the  computations  of  the  theoretically  derived  curves,  care  was  taken  to 
calculate  the  scattering/reflectivity  coefficient,  p,  for  each  plot,  in  accord¬ 
ance  with  theory.  To  this  end,  the  formula 

.  a  a  47tA/i  sin  \)/ 

Pav  =  e  ;  A<)>  =  - ^  -T 

was  applied,  where 
<)>  is  phase  lag  angle, 

A<j)  is  phase  difference  between  two  rays  reflecting  from  the  surface, 

A h  is  standard  deviation  of  normal  distribution  of  hill  heights, 

\jr  is  the  grazing  angle,  and 
X  is  wave  length. 
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5.4  Comparison  of  Theoretically  Derived  Power  Density  Loss 
for  Single  Diffraction  with  Data  for  Diffraction  over  Two 
Peaks 

In  appendix  G,  plot  G-l  shows  the  theoretically  derived  curve  of  power 
density  loss  due  to  diffraction  over  one  peak  (using  the  four-ray-path- 
theory  computer  program  describing  diffraction  over  a  single  edge,  with 
reflections  on  electrically  interactive  terrain  [3,7])  compared  to  field  meas¬ 
urement  data  for  signal  transmission  over  two  peaks. 

5.5  Plots  for  Doubly  Diffracted  Horizontally  Polarized  Micro- 
waves  Based  on  Theory  and  the  Application  of  Program 
DOUBLEPLOT 

It  should  be  noted  that  appendix  D  shows  plots  of  power  density  at  target 
heights  downrange  for  doubly  diffracted  horizontally  polarized 
microwaves. 

The  computer  program  DOUBLE_PLOT,  listed  in  appendix  C,  sets  up  the 
numerical  power  density  values  (as  the  abscissa)  versus  target  heights  (as 
the  ordinate)  for  each  plot  in  appendices  D  through  G. 
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Appendix  A.  —  Picturing  Diffractive  Ray-Path  Geometries  for 
Two  Knife  Edges 


The  following  diagrams  represent  unique  ray  paths  associated  with  each 
of  their  normalized  electric  amplitudes  Ajk,  where  j  =  2,...,  8,  and  k=  1,2. 

Note  that  in  this  appendix,  each  pair  of  ray  path  diagrams  (designated  An 
and  An)  shown  in  figure  A- 1  describes  the  diffraction  that  occurs  when  a 
ray  passes  over  two  knife  edges.  The  diagrams  of  the  figure  give  insight 
first  to  the  diffractive  influence  caused  by  the  passage  of  the  ray  over  peak 
F i  and,  second,  to  the  diffractive  influence  caused  by  passage  of  the  ray 
over  peak  Fz.  To  each  diagram  (An  or  Aa\  i  =  2, ....  8)  there  corresponds 
a  version  of  the  algebraic  expression  for  computing  the  total  indirect  ray 
path  amplitude  A,iAi2.  The  notation  designating  the  diagrams  of  the  fig¬ 
ure  and  the  notation  designating  the  respective  amplitudes  are  the  same. 
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APPENDIX  A 

Diffraction  over  F, 

Diffraction  over  F2 

Figure  A-l.  Indirect  ray  paths  and  their  corresponding  normalized  electric  amplitudes  (note:  path  1  Is 
the  direct  ray). 
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Appendix  B.  —  Computer  Program  DOUBLEDGE 
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APPENDIX  B 


Listed  below  are  the  input  parameters  of  the  computer  program  DOUBLEDGE: 


Geometric  parameters* 


Z1  is  antenna  height 

Z2S  is  minimal  receiver  height 

Z2E  is  maximal  hypothetical  receiver  height  of  interest 

DZ2  is  step  in  meters  to  next  height  level  of  receiver 

D1  is  distance  from  transmitter  to  first  knife  edge 

D2  is  distance  from  first  knife  edge  to  second  one 

D3  is  distance  from  second  knife  edge  to  receiver 

FI  is  height  of  first  knife  edge 

F2  is  height  of  second  knife  edge 


Electric  and  terrain  parameters 


EPSLN1 

EPSLN2 

EPSLN3 

51 

52 

53 

REFL1 

REFL2 

REFL3 

RLAMDA 

POLR 


is  £i,  which  is  the  relative  dielectric  constant**  over  D1 
is  £2,  which  is  the  relative  dielectric  constant*  over  D2 
is  £3,  which  is  the  relative  dielectric  constant**  over  D3 
is  conductivity  of  the  groundt  over  D1  in  mhos/meters 
is  conductivity  of  the  groundt  over  D2  in  mhos/meters 
is  conductivity  of  the  ground!  over  D3  in  mhos/meters 
is  the  roughness*  factor  over  D1 
is  the  roughness^  factor  over  D2 
is  the  oughnesst  factor  over  D3 
is  the  wavelength  in  meters 
is  the  polarization  of  the  electric  field  E  and  may  be 
either  horizontal  “H”  or  vertical  “V” 


AEOLOTROPIC  is  a  parameter  that  distinguishes  between  a  beamed  radar  signal  and  a 

broadcast/isotropic  lower  frequency  signal.  It  is  therefore  input  as  either 
“ANISOTROPIC”  or  as  “ISOTROPIC,”  with  an  appropriate  change  in  the 
wavelength  RLAMDA 


*All  heights  and  distances  are  in  meters. 

**This  is  also  called  permittivity. 

tThe  unit  mho  is  a  reciprocal  ohm.  One  ohm  is  the  same  as  one  siemen  in  SI  units. 
tThis  satisfies  the  Rayleigh  criterion. 
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APPENDIX  B 


C  PROGRAM  DOUHLEDGE 

r  THIS  MODEL  CALCULATES  THE  PATTERN  PROPAGATION  FACTOR  F  FOP  RADIO 

C  PROPAGATION  OVER  FLAT  TERRAIN  ON  WHICH  2  KNIFE-EDGE  OBSTRUCTIONS 

r  LIE  PERPENDICULAR  TO  THE  DIRECTION  OF  PROPAGATION. 

r  input: 

C  GEOMETRICAL  DISTANCES  IN  MFTERS 

C  WAVELENGTH  IN  METERS 

C  POLARIZATION:  H  OR  V 

C  GROUND  PROPERTIES  TO  THE  LEFT  AND  RIGHT  OF  EACH  MASK 

r  OUTPUT: 

C  -  TABLE  CONSISTING  CF: 

RADIATING  SOURCE  ANTENNA  HEIGHT  Z1  (IN  METERS) 

TARGET  HEIGHT  Z2  (METFPS) 

C  GRAZING  ANGLES  AT  THE  TRA'ISMITTFR  AND  THE  TARGET 

C  VARIABLES  THAT  START  WITH  THE  LETTER  C  ARE  COMPLEX 

CHARACTER*!  POLARC2  )  ,  PCLR 
CHARACTER* 1 1  BE AM ( 2  )  ,  AEOLOTROPIC 
COMPLEX  CT1  .CT2.CT3 .CF.CTH.CTV 
COMPLEX  CE11 ,  C£  1 2,CE21 , CE  3  1 ,  CF4  1 
COMFLEX  CE52.CE62.CE7l ,CE 72. CE 81 ,CE82 
COMPLEX  CE32 ,  CE42.  CE22,  CE51,  CE61 

REAL *8  R 1 ,R2.R3.R4,R5,R6,R7,R8,A.B, VTEPM.DD.D 

P  E A  L  *  8  VC11.VC12.VC21.VC31.VC41 

PEA  L  *8  VC52.VC62.VC71.VC72.VC81.VC82 

PE  A  L  *  8  VS11.VS12.VS21.VS31.VS41 

REAL*8  VS52.VS62,VS71,VS72.VSR1,VS82 

REAL*8  VC32.VC42.VC22.VC51 , VC61 

RE AL  *  8  VS32.VS42.VS22.VS51.VS61 

PEAL *4  PL A“D A 

DATA  POLAR/'H' , 'V'/ 

DATA  BEAM/ 'AN  ISOTROPIC', 'ISOTROPIC'/ 

C  OPEN  INPUT  AND  OUTPUT  FILES,  AND  DEFINE  SOME  OUANTITIF.S. 

OPEN  ( 1  ,  NAME  =  'DOURLEDGE.DAT',  TYPF  =  'OLD') 

OPE  N  (  2  ,  NAME  =  'DOUBLECGE .OUT' .  TYPE  =  'NEW') 

RE AD ( 1 , * )  LCASES 

PI  =  3.14159 
RAC  =  57.29578 
SORT2  =  SORT (2.0) 

AA  =  -82914.54 
PB  *  126.3854 
CC  =  48.07834 

RCS  =  1.  !  RADAR  CROSS- SFCTION 

DO  300  ICASE  =  1,  LCASES 
C  READ  GEOMETRIC  PARAMETERS. 

C  THREE  RECORDS  CONSTITUTE  A  SUP-FILE 
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r  pur  geometric  parameters 

C  Zl  IS  ANTENNA  HEIGHT 

C  Z2S  IS  minimal  RECEIVER  HEIGHT 

C  Z2E  IS  maximal  hypothetical  receiver  hftght  of  interest 

C  DZ2  IS  STEP  IN  METERS  TO  NEXT  HEIGHT  LEVEL  OF  RECEIVER 

C  D 1  IS  DISTANCE  FROM  TRANSMITTFP  TO  FIRST  KNIFE  EDGE 

C  D2  IS  DISTANCE  FROM  FIRST  KNIFE  EDGF  TO  THE  SECOND  ONE 

C  D 3  IS  DISTANCE  FROM  THE  SECOND  KNIFE  EDGE  TC  THE  RECEIVER 

C  FI  IS  HEIGHT  OF  FIRST  KNIFE  EDGE 

r  F2  IS  HEIGHT  OF  SECOND  KNIFE  EDGE 

READCl,*,  END  =  16)  Zl,  22 S,  Z2E,  DZ2 ,  Dl,  D2.  D3,  FI,  F2 

READ  IN  DIELECTRIC  CONSTANTS  AND  FRESNEL  FACTORS  AND  WAVELENGTH 
EPSLN1  IS  EPSILON  (PERMITIV1T Y )  -  RELATIVE  DIELECTRIC 
CONSTANT  OVER  Dl 

,  EPSLN2  IS  EPSILON  (PEPMITIVITY)  -  RELATIVE  DIELECTRIC 
CONSTANT  OVER  02 

EPSLN3  IS  EPSILON  ( FERMITIVITY )  -  RELATIVE  DIELECTRIC 
CONSTANT  CVEP  03 

51  IS  CONDUCTIVITY  OF  GROUND  OVER  Dt  IN  MHOS/METER 

52  IS  CONDUCTIVITY  OF  GROUND  OVER  D2  IN  MHOS/METER 

53  IS  CONDUCTIVITY  OF  GROUND  OVEP  D3  IN  MHOS/METER 

PEFL1  IS  ROUGHNESS  FACTOR  OVER  Dl 

REFL2  IS  ROUGHNESS  FACTOR  OVER  D? 

REFL3  IS  ROUGHNESS  FACTOR  OVER  D3 

PLAMDA  IS  WAVELENGTH  IN  METERS 

READ (  1 , * )  EPSLN1,  EPSLN2 ,  EPSLN3,  SI,  S2.  S3,  RFFL1, 

1  REFL2,  REFL3,  RLAMDA 

POLP  IS  POLARIZATION,  AND  MAY  EITHER  PF  VERTICAL  "V- 
OR  HORIZONTAL  *H". 

READ  IN  POLARIZATION 
READ  C 1 , 3  )  POLR 

ANISOTROPIC/AEOLOTRCPIC  MEANS  THAT  THE  MICROWAVES  ARE 
PF.AMED  S,  THEREFORE  NCT  BROADCAST  (I.E..  IS  TSCTROPTC) 

READ  Cl, 4)  AEOLOTROPIC 

TFCPOLR.EO.POLARCl  )  )  IPCL  =  1 

IFCP0LR.EQ.P0LARC2) )  IFCL  =  2 

IF (AEOLOTROPIC. EO. BEAM ( 1 )  )  ITPOPIC  =  1 
IF(AE0L0TR0PIC.EQ.BEAM(2) )  TTROPIC  =  2 

INITIALIZE  SOME  VALUES  BEFORE  THE  MAIN  PROCESSING  LOOP:  Z7  TS 
THE  HEIGHT  OF  RFCEIVER  (m),  LAMDA  IS  THE  WAVELENGTH. 

72  =  Z2S  -07.2 

r  FCHO  CHECK  THE  INPUTS  AND  PREPARE  THE  OUTPUT  TABLE. 

WRITEC2, 10)  Zl,  01,  D2,  03,  FI,  F2,  PLAMOA,  POLR,  EPSLNt, 

1  EPSLN?  ,  FPSLN  3 ,  SI,  S2 ,  S3,  REFL1,  REFL2,  REFL  3 

C  CALCULATE  THE  TWO  CLEARANCE  PARAMETERS 

C  FIRST  PARAMETER 

HOI  =  RLAMDA  »  Dl  *  02 


C 

C 

C 

C 

r 

r 

c 

r 

C 

C 

r 

r 

c 

r 
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ho2  =  m  *02 

HOO 1  =  SQRT(H01 /H02) 

r  SECOND  parameter 

H03  =  RLAMDA  *  D2  *  D3 

H04  =  02  ♦03 

H002  =  SQRT(H03/HQ4 ) 

r  MAIN  PROCESSING  LOOP:  LCOP  OVER  TARGET  HEIGHT  FROM  INITIAL 

C  TARGET  HEIGHT  TO  FINAL  TARGET  HEIGHT  BY  THE  INCREMENT  DZ2 

100  Z2  =  Z2  +DZ2 

IF ( Z2  .GT.  Z2E)  GOTO  1000 

C  COMPUTE  THE  GRAZING  ANGLES  FOR  CHANGES  OF  TARGET  HEIGHT 

TTHETA2  =  (Z1  +Z2)/(D1  +02  +  D3M.0D  +  00) 

TTHETA1  s  (Z2  -Z1)/(D1  +  D2  +  D3M.0D  +  00) 

PHI  s  TTHETA1  +TTHETA2 

IF( ITROPIC  .EO.  1)  GO  TC  102 
GRdB  =  1. 

I F ( ITROPIC  .EO.  2)  GO  TC  101 
102  CONTINUE 

G  =  AA*PHI*PHI  +BB+PHI  +  CC 

SUPER  =  (G  -CC)/20. 

GRdB  =  10. *  MSUPER 

C  CALCULATE  Rl,  R2,  R3.  AND  R4 

101  CONTINUE 

Rl  =  SORT ( ( Z2  -Zl)**2  >(01  +  D2  +D3)**2) 

R2  =  SORT ( C Z 1  +Z2)**2  +  (D1  +  P2  +D3)**2) 

R3  =  R? 

P4  =  Rl 
R5  =  R2 
R6  =  Rl 
P7  =  Rl 
PR  =  Rl 

C  CALCULATION  OF  DIFFRACT  ION  CLEARANCE  -  DELTA  PARAMETERS 


CALCULATE  Vll, 

V  1  2  , 

V21  , 

V  3 1  . 

V  4  1  , 

V  52, 

V62  , 

V71  , 

V  72, 

V  8  1  , 

V  8  2  , 

V61  , 

ALSO  CALCULATE 

V  32 , 

V42  , 

V22  , 

V  4  2  , 

V51 

0  THE  END  CF  THIS 

C  GROUPING  OF  INSTRUCTIONS 

D  =  D1  >02*1 .00+00 
00  =  D  +03 

C  CLEARANCE  DELTA  IS  (Z1  +VTERM  -FI)  &  DFLTA-ZERC  IS  H001 

VTF.RM  =  (Z2  -Z  1  )  *D  1  /DO 
Vll  =  S0RT2  *  ( Z 1  +VTERM  -FD/H001 
1  INTERPRETATION.  DIRECT  CASE. 

C  CLEARANCE  DELTA  IS  (Z1  +VTERM  -F2)  4  DELTA-ZERO  IS  H002 

VTEPM  =  (Z2  -Z1)*0/DD 

V 1 2  =  SQRT2  *  (Z1  +VTERM  -F2)/H002 

!  INTERPRETATION,  DIRECT  CASE. 

C  CLEARANCE  DELTA  IS  (-Z1  +VTERM  -El)  t  PELTA-ZERC  IS  H001 

VTEPM  =  (Z1  +Z2)*D1/DD 
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V?1  =  SQWT2  *  (-Z1  +  VT  FP*  -FD/HOOl 

r  C  LFARANCE  DELTA  IS  C  Z 1  +VTERM  -El)  5,  DELTA-ZFRC  IS  HCOl 

VTERM  =  (-Z1  -7.2)*01/DC 
V 3 1  =  SQRT2  *  CZ1  ♦VTERM  -FD/HOOl  . 

C  CLEARANCE  DELTA  IS  (-Z1  +VTFRM  -El)  £.  DFLTA-ZEPO  IS  HOOl 

VTEPM  =  (21  -Z2)*D1/DD 
V41  =  SORT2  *  (-Z1  +VTEPM  -FD/HOOl 

C  CLEARANCE  DELTA  IS  (Z1  ♦VTERM  -E2)  &  DELTA-ZERO  IS  H002 

VTERM  =  (-Z1  -Z2)*D/DD 
V52  =  SQRT2  *  (Z1  +VTERM  -F2)/H002 

r  CLEARANCE  DELTA  IS  (-Z1  +VTFRM  -E2)  5,  DELTA-ZFRO  IS  H002 

'  VTERM  a  (Z1  -Z2)*D/DD 

V62  a  SQRT2  *  (-Z1  eVTERM  -F2)/H002 

C  CLEARANCE  DELTA  IS  (-Z1  +VTERM  -FD  i  DFLTA-ZEPC  IS  HOOl 

VTEPM  a  (Z1  -Z2)*D1/0D 
V71  a  SORT2  *  (-Z1  ♦  V  T  E  R  M  -FD/HOOl 

C  CLEARANCE  DELTA  IS  (-Z1  +  VTERM  -F2)  t  DFLTA-ZERO  IS  H002 

VTEPM  a  (Z1  -Z2)*D/DD 

V72  a  SORT2  *  (-Z1  +VTERM  -F2)/H002 

r  CLEARANCE  DELTA  IS  (-Z1  +VTEPM  -FD  6  DFLTA-ZERO  IS  HOOl 

VTERM  a  (-7.2  ♦  ZlDDl/DC 
V81  a  SOPT2  *  (-Z1  +VTERM  -FD/HOOl 

C  CLEARANCE  DELTA  IS  (-Z1  ♦ VTERM  -F2)  &  DFLTA-ZFRO  IS  HC02 

VTERM  =  (-Z2  +Z1  )*D/DD 
V82  =  SQRT2  *  (-Z1  ♦  V  T  E  R M  -F2)/H007 

VTERM  =  (Z1  ♦Z?)*P/nD 

V 32  =  SORT2  *  (-Z1  +VTERM  -F2)/H002 

VTERM  =  (Z1  +Z2)*D /on 

V42  =  SORT2  ♦  (-Z1  ♦VTEPM  -F?)/H002 

VTERM  =  (Z1  ♦Z2)*D/nD 

V22  =  SORT2  *  (-Z1  ♦VTEPM  -F?)/H002 

VTEPM  =  (-71  - Z?  )  * D  1  /  DC 

V51  =  SOPT2  ♦  (Zi  +VTERM  -Fl)/Hn01 

VTERM  a  (-Z1  -Z2)*D1/DC 

V61  =  SQRT2  *  (Zl  ♦  VT  E  P  M  -FD/HOOl 

0  TEST  FOR  CHOICE  OF  CIRFOT  PATH  PAY 

ALPHA1  =  AT  AN  (  (FI  -ZD/D1  ) 

ALPHA2  =  AT  AN  (  C  F  2  -ZD/D  ) 

C  TF(  (FI.  LT.  F 2)  .AND.  (ALFHAt  .LT.  ALPHA2  )  )  THEN 

C  GO  TO  52 

c  FLSE 

c  GO  TO  50 

C  FNP  IF 

C  I F (  (FI  .LT.  F2)  .AND.  ( A  L  PH  A 1  . GE .  ALPHA2)  )  THEN 
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c  GCTO  SO 

C  Ft  SF 

c  GC  TO  52 

C  FND  IF 


c 

IF  ( 

(FI  .GE.  F2 )  .AND.  (ALF  <A1 

.GE.  ALPHA2) 

) 

THEN 

c 

GO  TO  50 

c 

FLSE 

c 

GC  TO  52 

c 

end 

IF 

IF  ( 

(FI  .LT.  F2)  .AND.  ( ALPHA1 

.LT.  ALPHA2 ) 

) 

GOTO 

52 

IF  ( 

(FI  .LT.  F2 )  .AND.  ( ALPHA  1 

.GE.  ALPHA2) 

) 

GOTO 

50 

IFC 

(FI  .LT.  F2 )  .AND.  ( ALPHA1 

.GE.  ALPHA2 ) 

) 

GOTO 

50 

C  CALCULATE  All  USING  VII  CORRESP.  TO  FI 

so'  call  ncscvcn,  VS11,  Vll) 

All  =  1./SQRT2  *  DSGPT ( ( VC  1  1  +0.5)**2  +  (VS  11  +0.5)**2) 

TFCVCtl  +0.5  .GE.  0.)  THEN 

B  =  DATAN(  (VS11  +  0.5)/(VCU  +0.5)  )  -PI/4. 

ELSE 

B  =  PI  +OATANC  (VS11  +0.5)/(VCll  +0.5)  )  -PI/4. 

END  IF 

A 1  =  All  *  DSIN(e) 

B 1  =  All  *  occsce) 

CE11  =  CMPLX(ei.Al) 

CE1  =  CE1 1 
GOTO  53 

C  CALCULATE  A12  USING  V12  CORRFSP.  TO  F2 

5?  CALL  DCS( VC12  #  VS12,  V12) 

A  1 2  =  1./S0RT2  ♦  OSORT ( { VC  1 2  +0. 5)4*2  +(VS12  +0. 5)**2) 

I F ( VC  1  2  +0.5  .GE.  0.)  THEN 

B  =  DAT AN ( (VS12  +0.5)/(VC12  +0.5))  -PI/4. 

ELSE 

B  =  PI  +DATANC (VS1?  +0.5)/CVCl?  +0.5))  -PI/4. 

END  IF 

A 1  =  A 1 2  *  DSIN(B) 

B 1  =  A 1 2  *  DCCS(e) 

CE12  =  CMP[,X(ei,Al) 

CE1  =  CE 1 2 

r  COMPUTATION  OF  GRAZING  ANGLE  PSI1  OVEP  D1  , 

r  PS  12  OVER  02,  AND  PSI3  CVEP  03. 

C  GRAZING  ANGLE  FOR  PATH  21 

S3  IFCV21  .GE.  0.)  GOTO  200 

C  PAY  INTERSECTS  THE  BASF  FI 

PSI1  =  D  A  T  A  N (  C  Z 1  +F1)/D1*1 .OD+OO) 

GOTO  201 

r  PAY  IS  CLEAR  OF  THE  MASK  FI 

200  CONTINUE 

PSI1  =  DAT AN (  (Z2  +Z1)/CD) 
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201  CONTINUE 

CALL  THE  COMPLEX  REFLECTION  SUBROUTINE  FOR  THE 
TERRAIN  REGION  D1 

CALL  FRESNLCEPSLNl .RLAMCA.Sl ,PSI1 ,CTH,CTV) 

C  MULTIPLY  BY  FACTORS 

CTH  =  CTH  *  REFL1 
CTV  =  CTV  *  REEL  I 

C  SET  CT1  OVER  D! 

IE ( I  POL  .EQ.  1)  CT1  =  CTH 
IFUPOL.  EO.  2)  CT1  =  CTV 

C  888  SHOWS  NO  REFLECTION  OVER  D1 

IF(EPSLN1  .EQ.  888.  .AMD.  SI  .EO.  888.)  CT1  =  CMPLX f 0 . 0 . 0 . 0 ) 

C  CALCULATE  A21  USING  V21  CORRESP.  TO  El 

CALL  DCSIVC21 ,  VS21,  V21) 

A21  =  1./SQRT2  »  DSGRT ( ( VC 2 1  +0.5)**2  +IVS21  +0.5)**2) 

C  GRAZING  ANGLE  FOR  PATH  22 

IFCV22  .GF.  0.)  GOTO  220 

C  RAY  INTERSECTS  THE  MASK  F2 

PSI1  =  DATA  N (  (Z1  +F2)/C) 

GOTO  221 

C  PAY  IS  CLEAR  OF  THE  MASK  F2 

220  CONTINUE 

PSI1  =  D AT  A  N (  (Z1  +Z2)/C0) 

221  CONTINUE 

CALL  F PFSNL(EPSLN1 , RL A M C A . S 1 , PS I  1 , CTH,  CTV) 

CTH  =  CTH  *  RFFL1 

CTV  =  CTV  *  REFL1 

IF ( I  POL  .EO.  1)  CT1  =  CTH 

I F  C I  POL  .EQ.  2)  CT1  =  CTV 

I F ( E PSLN 1  . EQ .  888.  .AND.  SI  .EQ.  888.)  CT1  =  CMPLX ( 0 . 0 , 0 . 0 ) 
C  CALCULATE  CE27  USING  V72  CORRESP.  TO  E2 

CALL  DCS ( VC22 ,  VS 22,  V22) 

A  22  =  1./SQRT2  *  DSQRT  C ( VC2  2  +0.5)**2  +CVS22  +0.5)**2) 

IF ( VC??  +0.5  .GE.  0.)  THEN 

B  =  DATAN((VS?2  +0 , 5 ) / ( VC 2 ?  +0.5))  -PI/4. 

ELSE 

e  =  PI  +DATAN((VS22  +0.5)/(VC22  +0.5))  -PI/4. 

END  IF 

A 1  =  A?2  *  A  2 1  *  QSIN(P) 

B 1  =  A  22  *  A  2 1  *  DCOS(P) 

CF22  =  CMPLX  (PI ,  A 1  )  *  CT1  *  GRrJB 

0  GRAZING  ANGLE  FOP  PATH  31 

I F ( V  3 1  .GE.  0.)  GOTO  310 
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r  PAY  INTERSECTS  THE  MASK  El 

psn  =  n  a  t  a  n  (  cz?  +fi)/(02  +m*i.on+oo)  ) 

GOTO  301 

C  RAY  IS  CLEAR  CF  THE  MASK  FI 

310  CONTINUE 

PSI1  =  DATA  N (  (71  +Z2)/DD) 

301  CONTINUE 

CALL  FRESNL(EPSLN1 ,  RL AMD A . S 1 , PS  1 1 , CTH . CTV ) 

CTH  =  CTH  *  REFL1 

CTV  =  CTV  *  REFL1 

I F ( I  POL  .EG.  1)  CT1  =  CTH 
I F ( I POL  .EG.  2)  CT1  =  CTV 

IF (EPSLN1  .EG.  888.  .AND.  SI  .EQ.  888.)  CT1  =  CMPLX ( 0 . 0 . 0 . 0  ) 
C  '  CALCULATE  A31  USING  V31  CORRESP.  TO  FI 

CALL  DCS(VC31,  VS  3  1  ,  V31) 

A  3 1  =  1./S0RT2  *  DSQRT ( ( VC3 1  +0.5)**2  +CVS31  +0.5)**2) 

r  GRAZING  ANGLE  FOR  PATH  32 

I F ( V 32  .GE.  0.)  GOTO  320 

C  RAY  INTERSECTS  THE  MASK  F2 

PSI1  =  D AT  A  N (  (Z1  +F2)/D) 

GOTO  321 

C  RAY  IS  CLEAR  OF  THE  MASK  F2 

370  CONTINUF 

PSI1  =  DATANC  (Z1  +Z2 }/DP) 

321  CONTINUE 

CALL  FRESNUEPSLNl  ,RLAVDA,S1  ,PSI1  , CTH, CTV) 

CTH  =  CTH  ♦  REFL1 

CTV  =  CTV  *  REFL1 

I F ( I  POL  .EG.  1)  CT1  =  CTH 
IFCIPOL.  EQ.  2)  CT1  =  CTV 

IFCEPSLNl  .EQ.  888.  .AND.  SI  .EG.  6  8  e  .  )  CT1  =  CMPLX ( 0 .0 , 0 . 0  ) 
C  CALCULATE  CE32  USING  V32  CORPESP.  TO  F7 

CALL  DCS ( VC 32 ,  VS32,  V37) 

A  32  =  1./SQRT2  *  DSOPT((VC32  +0.5)**2  +CVS32  +0.5)**2) 

I E ( VC  32  +0.5  .GE.  0.)  THFN 

B  =  DATAN( (VS32  +0.5)/(VC32  +0.5))  -PI/4. 

ELSE 

R  =  PI  +DA7ANC  (VS32  +0.5)/(VC32  +0.5))  -PI/4. 

END  IF 

A 1  =  A32  *  A  3 1  ♦  DSIN(B) 

B 1  =  A32  *  A3 1  *  DCOS(B) 

CE 32  =  CMPLX(P1,A1)  »  CT2  *  CRrlR 

C  GRAZING  ANGLE  FOR  PATH  41 

IF ( V 4 1  .GE.  0.  )  GOTO  400 

C  RAY  INTERSECTS  THE  MASK  FI 
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PSU  =  D  A  T  A  N  (  (FI  4Z1)/  HI  *  l.OD  +  OO  ) 

PS1 2  =  DATAN (  ( Z2  +F1)/  0  ) 

CALL  FRESNL ( EPSLN 1 .RLAMOA.Sl , PSI1 .CTH.CTV) 

CTH  =  CTH  *  REFLI 

CTV  =  CTV  *  REFL1 

IF ( IPOL  .EO.  1)  CT1  =  CTH 

IF ( IPOL  .EO.  2)  CT1  =  CTV 

IF (EPSLN 1  .EQ.  888.  .ANC.  SI  . EO.  888.)  CT1  =  CMPLX ( 0 . 0 , 0 . 0 ) 

C  CALL  THE  COMPLEX  COEFFICIENT  OVER  D2 

CALL  FRESNL(EPSLN2,RIAMDA,S2.PSI2,CTH,CTV) 

C  MULTIPLY  BY  FACTORS 

CTH  =  CTH  *  REFL2 

,  CTV  s  CTV  *  REFL2 

C  SET  CT2  OVER  02 

I F ( IPOL  .EO.  1)  CT2  =  CTH 
IF ( IPOL  .EO.  2)  CT2  =  CTV 

r  888  SHOWS  0  REFLECTION  OVER  D2 

IF ( EPSLN2  .EO.  888  .  .ANC.  S2  .EO.  888  .)  CT2  =  CMPLX ( 0 . 0 , 0 . 0 ) 

GOTO  401 

C  RAY  IS  CLEAR  OF  THE  MASK  FI 

400  CONTINUE 

PSI1  =  DATANf  C Z 1  -Z2 ) /  DC  ) 

CALL  FRFSNL( EPSLN 1 .RLAMOA.Sl ,PSI1 , CTH, CTV) 

CTH  =  CTH  *  REFLI 
CTV  s  CTV  *  REFLI 
I F ( IPOL  .EO.  1)  CT1  =  CTH 
I F ( IPOL  .EQ.  2)  CT1  =  CTV 

I F ( EPSLN 1  , EO .  888  .  .ANC.  SI  .EQ.  888  .)  CT1  =  CUPLX ( 0 . 0 , 0 . 0 ) 

401  CONTINUE 

C  CALCULATE  A41  USING  V41  CCRRESP.  TO  FI 

CALL  DCS ( VC4 1  ,  VS41,  V41) 

A 4 1  s  1./S0RT2  *  DSQFT ( ( VC 4 1  ♦  0.5)**2  ♦  (VS41  4  0.5)**2) 

r  GRAZING  ANGLE  FOR  PATH  42 

I F ( V 42  ,GE.  0.  )  GOTO  420 

r  RAY  INTERSECTS  THE  MASK  F2 

PSI1  =  DATAN (  (Z1  +F1)  /  01  *  1.0D400  ) 

C  PSI 2  =  DATAN (  (Z1  4E2)/C) 

CALL  FRESNL ( EPSLN 1 , PL  AMD A, SI ,PSI1 .CTH.CTV) 

CTH  s  CTH  *  REFLI 
CTV  =  CTV  *  REFLI 
IF( IPOL  .EQ.  1)  CT1  =  CTH 
TF( IPOL  .EO.  2)  CT1  =  CTV 

I F ( EPSLN 1  .EO.  888  .  .AND.  SI  .EO.  888  .)  CTl  a  C M PL X ( 0 . 0 , 0 . 0 ) 

CALL  FRESNL (EPSLN2.PL AMD A ,S2, PSI 2, CTH, CTV) 

CTH  s  CTH  *  REFL2 
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CTV  =  CTV  *  PEEL2 
1 F ( I  POL  .EO.  n  CT2  =  CTH 
I F ( I  POL  .EO.  2)  CT?  =  CTV 

IF (EPSLN2  .EO.  888.  .AND.  S2  .EO.  888  .)  CT  2  =  CM PLX ( 0 . 0 , 0 . 0  ) 
GOTO  421 

C  PAT  IS  CLEAR  OF  THE  MASK  F2 

420  CONTINUF 

PSI1  =  DATAN (  (Z1  +Z2)/DD) 

CALL  FRESNL(EPSLN1 ,RLAMDA,S1 .PSI1 , CTH, CTV) 

CTH  =  CTH  *  REFL1 
CTV  =  CTV  *  REEL  1 
I F ( I  POL  .EO.  1)  CT 1  *  CTH 
I F ( I  POL  .EO.  2)  CT1  =  CTV 

IF ( EPSLN 1  .EQ.  888.  .AND.  SI  .EO.  888  .)  CT1  =  CMPLX (0.0 , 0 .0) 

/ 

421  CONTINUE 

C  CALCULATE  CE42  USING  V42  CORRESP.  TO  F2 

CALL  DCS ( VC42 ,  VS42,  V42) 

A42  =  1  .  /SORT  2  *  DS0RTUVC42  ♦0.5)**2  +CVS42  +0.5)**2) 

IF( VC4  2  *0.5  .GE.  0.)  THEN 

H  =  OATAN ( ( VS42  +0 . 5 ) / ( VC42  +0.5))  -PI/4. 

ELSF 

8  =  PI  ♦DATAN((VS42  +0.5)/(VC42  +0.5))  -PI/4. 

END  IF 

ftl  =  A 4 2  *  A 4 1  *  DSIN(R) 

PI  r  A4?  *  A 4 1  *  DCOS(P) 

CE42  =  CMPLX(Bl.Al)  *  CT1  »  CT2  *  GRdP 

C  GRAZING  ANGLE  FOR  PATH  51 

I  F ( V  5 1  .GE.  0.  )  GOTO  590 

C  RAY  INTERSECTS  THE  MASK  FI 

PSI3  =  DATANC  ( Z  2  +F2  )  /  n3»1.0C  +  00) 

C  CALL  THE  COMPLEX  COEFFICIENT  OVER  D3 

CALL  FRESNLCEPSLN3 ,RLAMOA ,S3,PSI 3 ,CTH ,CTV) 

C  MULTIPLY  BY  FACTORS 

CTH  =  CTH  *  REFL3 
CTV  =  CTV  *  RFFL3 

C  SET  CT 3  OVER  D3 

I  F ( I  POL  .EO.  1)  CT  3  =  CTH 

I F (  I  POL  .FQ.  2)  CT  3  =  CTV 

C  IF(EPSLN3  .EO.  888.  .AND.  S3  .EQ.  888  .)  CT  3  =  CM  PL  X ( 0 . 0 , 0 . 0  ) 

GOTO  501 

C  RAY  CLEAR  OF  THE  MASK  FI 

590  CONTINUE 

PSI1  =  DATANC  (Z1  +Z2)/  DD) 

CALL  ERESNL( EPSLN 1 ,PLAMDA,51,PSI1,CTH,CTV) 

CTH  =  CTH  *  REFL1 
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CTV  =  CTV  *  R  E  F  t  >  1 
I F  ( I POL  .EO.  1)  CT1  =  CTH 
I F  ( I  POL.  .EG.  2)  CT1  =  CTV 

I F ( FPSLN 1  .EQ.  888.  .AND.  SI  .EG.  888.)  CT1  =  CMPLX ( 0 . 0 , 0 . 0 ) 
SOI  CONTINUE 

C  CALCULATE  A51  USING  V51  CORRESf*.  TO  FI 

CALL  DCS( VC51 ,  VS51,  V51) 

A51  =  1 . /SORT 2  *  OSQRT ( ( VC51  +0. 51**2  +(VSS1  +0.51**21 

C  GRAZING  ANGLE  FOR  PATH  52 

IF ( V 52  ,GE.  0.  1  GOTO  520 

C  RAY  INTERSECTS  THE  MASK  F2 

PSI 3  =  DATANC  CZ2  +F21/D3  *  l.OD+00  1 
CALL  FRESNL(EPSLN3,PLAMDA , S3, PSI 3 , CTH, CTV) 

'  CTH  =  CTH  *  REFL3 

CTV  =  CTV  *  REFL3 
I F ( I  POL  .EO.  1 1  CT  3  =  CTH 
IF  C I  POL  .EO.  2)  CT  3  =  CTV 

IF  CEPSLN  3  .EG.  888  .  .AND.  S3  .EO.  88  8.1  CT  3  =  CMPLX (0 . 0 , 0 . 0 ) 
GOTO  521 

RAY  IS  CLEAR  CF  THF  MASK  F2 

CONTINUE 

PSI1  =  DATAN (  (Z1  ♦ Z  2 1 /  DD  1 
CALL  FRESNLIEPSLNl , P L A MD A . S 1 , PS  1 1 , CTH, CTV ) 

CTH  =  CTH  *  REEL  1 
CTV  =  CTV  *  D EFL 1 
I F  ( I  POL  F  ).  11  CT 1  =  CTH 

I F  (  I  POt,  .tO.  21  CT1  =  CTV 

I F ( r  P  c :  .1  .EO.  808  .  .AND.  SI  .  EO  .  888.1  CT1  s  CMPLX ( 0 . 0 , 0 . 0 1 
CONTINUE 

CALCULATE  CE5 2  USING  V52  CORRESP.  TO  F2 

NOTE.  DITTO  AS  ABOVE  FOR  CASE  51.  R5  -R1  *  0 

CALL  DCS(VC52,  VS52,  V521 

A  5  2  =  1  . /SORT  2  *  DSOPT((VC52  +0.5)**?  +CVS52  +n.5)**?l 

IFCVC52  +  0.5  .GE.  0.1  THEN 

B  =  DATAN( ( VS52  +0.5)/(VC52  +0.511  -PI/4. 

ELSE 

e  a  PI  +DAT  AN ( ( VS52  +0.51/CVC52  +o.5)l  -°I/4. 

END  IF 

B  =  B  + ( 2 .  *  PI)/  PLAMDA  *  ( R5  -R 1  ) 

Al  =  A51  *  A52  *  DSIN(B) 

B 1  =  A52  *  A51  *  DCOSCB) 

CE52  =  CMPLX  (  B 1  .All  *  CT  3  *  GRrJB 

C  511  GRAZING  ANGLE  FOR  PATH  61 

I F (  V61  .GE.  0.  1  GOTO  6000 

C  RAY  INTERSECTS  THE  MASK  FI 

PS  1 2  =  DATAN (  (22  +F1 )  /  (D2  +D3*1.0D+00  )  1 
PS  1 3  =  DATANC  (7.2  +E2)  /  C3  *  l.OD+OO  1 

CALL  FRESNL(EPSLN2,RIAM0A ,S2 ,PSI2 ,CTH .CTV) 

CTH  s  CTH  *  REFL2 
CTV  «  CTV  *  REFL2 
IFCIPOL  .EQ.  1)  CT2  =  CTH 


r 

520 


521 
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IFCIPOL  .KG.  2)  CT  2  =  CTV 

IF (  E  P  5  L  N  2  .FO.  888.  .AMO.  $2  .  FQ .  8PP.)  CT2  =  C «PLX ( 0 . 0 , 0 . 0 ) 

CALL  FRESNLCEPSLN3,PLAMDA,S3,PSI 3,  CTH. CTV) 

CTH  =  CTH  *  REFL3 
CTV  =  CTV  *  REFL3 
IFCIPOL  .EO.  1)  CT3  =  CTH 
IFCIPOL  .EO.  2)  CT3  =  CTV 

IFCEPSLN3  .EO.  888.  .AND.  S3  .FQ.  888  .)  CT  3  =  CMPLX C 0 . 0 , 0 . 0  ) 
GOTO  601 

C  RAY  CLEAR  OF  THE  MASK  FI 

6000  CONTINUE 

PSI1  =  DATANC  CZ1  +Z21/CD  ) 

CALL  FRESNLCEPSLH1 .RLAMDA.Sl ,PSI1 ,  CTH,  CTV) 

'  CTH  =  CTH  *  REFL1 

CTV  =  CTV  *  REFL1 
IFCIPOL  .EC.  1)  CT 1  =  CTH 
IFCIPOL  .EO.  2)  CT 1  =  CTV 

IFCEPSLN1  .EO.  888  .  .AND.  SI  .EO.  88  8  .)  CTl  =  C MPLX ( 0 . 0 , 0 . 0  ) 
601  CONTINUE 

C  CALCULATE  A6l  USING  V61  CORRESP.  TO  FI 

CALL  DCSCVC61,  VS61,  V61) 

A 6 1  =  1./S0RT2  *  0SQRTCCVC61  -»0.5)**2  +CVS61  +0.5)**2) 

r  GRAZING  ANGLE  FOR  PATH  62 

IFCV62  .GE.  0.  )  GOTO  620 

r  HAY  INTERSECTS  THF  MASK  F2 

PS  1 2  =  DATANC  CZ1  +  F2)  /  D  ) 

PSI3  =  DATANC  CF2  +Z2)/C3  *  l.OD+OO  ) 

CALL  FRESNLCEPSLN2,PLAMDA,S2.PSI2,CTH,CTV) 

CTH  =  CTH  *  REFL2 
CTV  =  CTV  »  REFL2 
IFCIPOL  -EG.  1)  CT2  =  CTH 
I F ( I  POL  .EO.  2)  CT2  =  CTV 

IFCEPSLN2  .EO.  888  .  .AND.  S2  . EO .  8R8.)  CT2  =  CMPLX C 0 . 0 , 0 . 0  ) 

CALL  FRESNLCEPSLN3 , PL  AMD A  .S3.PSI 3 , CTH , CTV) 

CTH  =  CTH  *  REFL  3 
CTV  =  CTV  ♦  RFFL  3 
I F ( T  POL  .EQ.l  )  CT  3  =  CTH 
IFCIPOL  .EO.  2)  CT 3  =  CTV 

I F  C  FPSLN  3  .EO.  888.  .AND.  S3  .EQ.  R88.)  CT  3  =  CMPLX ( 0 . 0 , 0 . 0 ) 
GOTO  621 

C  RAY  CLEAR  OF  THF  MASK  F2 

620  CONTINUE 

PSI1  =  DATANC  CZ1  -Z2)/  DD  ) 

CALL  FRESNLCEPSLN1 .RLAMDA.Sl, PS  II , CTH, CTV) 

CTH  =  CTH  *  REFL 1 
CTV  =  .CTV  ♦  REFL  1 
IFCIPOL  .EO.  1)  CTl  =  CTH 
IFCIPOL  .EO.  2)  CTl  a  CTV 

IFCEPSLN1  .EO.  888.  .ANO.  51  .EO.  888.)  CTl  =  C MPLX C 0 . 0 , 0 . 0  ) 

621  CONTINUE 
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r 

r 


c 

* 

c 

r 

7000 

701 


C 

c 


r 

720 


CALCULATE  CEh 2  USING  V62  CCRRfSP.  TO  F2 

NOTE.  DITTO  AS  ABOVE  FCR  CASF  52.  R6  -R1  =  0 

CALL  DCS ( VC62 ,  VS 62,  Vt>2) 

A62  =  1./S0RT2  *  DSQRT((VC62  +0.5 )**2  +(VS62  +0.5)»*2) 

IF( VC62  +0.5  .GE.  0.)  THEN 

B  =  OATAN ( ( VS62  +0.5)/(VC62  +0.5)1  -PI/4. 

ELSE 

P  =  PI  +DATAN ( ( VS62  +0.5)/(VC62  +0.5))  -PI/4. 

END  IF 

Al  =  A  6 1  *  A62  *  DSIN(E) 

Pi  =  A62  *  A6 1  *  OCOS(B) 

CE62  =  CMPLX (  B 1  ,  A  1 )  *  C72  *  CT 3  ♦  GRdB 

GRAZING  ANGLE  FOR  PATH  71 
IF  ( V 71  .GE.  0.)  GOTO  7000 

PAY  INTERSECTS  THE  MASK  FI 

PSI1  =  DAT AN (  (FI  +Z11/D1  *  l.OD+OO  ) 

GOTO  701 

RAY  CLEAR  OF  THF  “ASK  FI 
CONTINUE 

PSI1  =  DAT  AN (  (Z1  -22)/  DD  ) 

CONTINUF 

CALL  FRESNLCEPSLNl , PL AMD A, SI, PS  II, CTH.CTV) 

CTH  =  CTH  *  REFL1 
CTV  =  CT V  *  REFL1 
I F ( I  POL  .EG.  1)  CT  1  =  CTH 
IF ( I POL  .EO.  2)  CT 1  =  CTV 

IF  l  FP5LN  1  .EQ.  989.  .  AND.  SI  .F.C.  888.)  CT1  s  CMPLX  (  0 . 0 , 0 . 0  ) 

CALCULATE  A7l  USING  V71  CORRESP.  TO  FI 
NOTE.  DITTO  AS  ABOVE  FOR  CASE  51.  R7  -R1  =  0 
CALL  DCS (VC71.VS71.V71) 

A  7 1  =  1./S0RT2  *  DS0PT((VC7)  +0.5)»*2  +(VS71  +0.5)**2) 

GRAZING  ANGLE  FOR  PATH  72 
IFCV72  .GE.  0.)  GOTO  720 

RAY  INTERSECTS  THE  “ASK  F2 

PS  1 3  =  DAT AN (  ( Z  2  +F2)  /  D3  *  l.OD  +  OO  ) 

CALL  FRESNL( EPSLN3.PL A “DA. S3, PS I  3, CTH.CTV) 

CTH  =  CTH  ♦  REFL3 

CTV  =  CTV  *  REFL3 

I F ( I  POL  .EO.  1)  CT  3  =  CTH 
I F ( I  POL  .EG.  2)  CT  3  =  CTV 

I F ( EPSLN  3  .EO.  888.  .AND.  S3  .EO.  888.)  CT 3  =  C«PLX ( 0 . 0 , 0 . 0  ) 
GOTO  721 

RAY  CLEAR  OF  THE  “ASK  F2 
CONTINUE 

PSI1  =  OATAN (  (71  -Z21/DD  ) 

CALL  FRESNL(EPSLN1 ,RL AMD A.Sl , PS II .CTH.CTV) 

CTH  s  CTH  *  REEL  1 

CTV  s  CTV  *  REFL1 

IFflPOL  .EG.  1)  CT  1  =  CTH 
I F ( I  POL  .EO.  2)  CT1  =  CTV 
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IFCEPSLNl  .EQ.  8  8  8.  .AMD.  SI  .FO.  8»8.)  CT1  =  C * P L X ( 0 . 0 , 0 . 0 ) 
721  CONTINUE 

CALCULATE  CE72  USING  V72  CORRESP.  TO  F2 
NOTE.  DITTO  AS  ABOVE  FOR  CASE  51.  R?  -R1  =  0 
CALL  DCS ( VC7  2  ,  VS72,  V72) 

A72  =  1 . /SORT 2  ♦  OSGRT C ( VC72  ♦0i5)**2  +CVS72  +0.5)**2) 

IFCVC72  +0.5  .GE.  0.)  THEN 

B  =  D ATAN  C ( VS7  2  +0.5)/(VC72  +0.5))  -PI/4. 

ELSE 

B  =  PI  +DAT  ANC (VS72  +0.5)/(VC72  +0.5))  -PI/4. 

END  IF 

A 1  =  A72  *  A71  *  DSIN(B) 

B 1  =  A72  *  A71  *  OCOSCB) 

'  CE72  =  CMPLXCei.Al)  *  CT1  *  CT 3  *  GRdB 

C  GRAZING  ANGLE  FOR  PATH  81 

I F  (  V  8  1  .  C.E  .  0.)  GOTO  800 

C  RAY  INTERSECTS  THE  MASK  FI 

PSI1  =  DATANC  (  21  +F1)/D1  »  l.OD+OO  ) 

PSI 2  =  DATANC  (FI  +Z2)/(02  +D3  *  l.OD+OO  )  ) 

CALL  FRESNLCEPSLNl .RtAMDA.Sl.PS II, CTH, CTV) 

CTH  =  CTH  *  REFL1 
C TV  =  CT  V  *  REFL1 
IF ( I  POL  .EO.  1)  CT 1  =  CTH 
I F ( I  POL  .EO.  2)  CT 1  =  CTV 

I  F ( FPSLN 1  .EO.  888.  .AND.  SI  .EO.  RRR.)  CT1  i  CMPLX ( 0 . 0 , 0 . 0  ) 

CALL  FRESNLCEPSLN2.RLAMnA.S2.PSI2.CTH.CTV) 

CTH  =  CTH  *  REFL2 
CTV  =  CTV  ♦  REFL2 
I F ( I  POL  .EO.  1 )  CT2  =  CTH 
I F ( I  POL  .FO.  2)  CT  2  =  CTV 

IF(FPSLN2  .EC.  888  .  .AND.  S2  .EO.  8RP.)  CT  2  =  C M P L X ( 0 . 0 , 0 . C  ) 
GOTC  801 

C  RAY  CLEAR  OF  THE  MASK  FI 

800  CONTINUE 

PSI1  =  DATANC  (F2  +Z2 ) /C  3  *  l.OD  +  OO  ) 

801  CONTINUE 

C  CALCULATE  A81  USING  V81  CORRESP.  TO  FI 

CALL  DCSCVC81,  VS8 1 ,  V81) 

API  =  1./SQRT2  ♦  DSQRTCCVC81  +0.5)**2  ♦  (  V  S  8  1  +0.5)*«2) 

C  GRAZING  ANGLE  FOR  PATH  82 

IF ( V 82  , GE .  0.)  GOTO  820 

C  RAY  INTERSECTS  THE  MASK  F2 

PSI  3  =  DATANC  ( F2  +Z2)/  03  *  l.OD  +  OO  ) 

CALL  FRESNLCEPSLN3  ,Pl  AMOA  ,'S  3  ,  PS  I  3  ,  CTH  ,  CT  V  ) 

CTH  s  CTH  *  REFL3 
CTV  s  CTV  *  REFL3 
I F ( I  POL  .EO.  1 )  CT  3  =  CTH 
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c 

820 


821 

C 


c 


r  215 

r 

216 


700 


I F ( I  POL  .EO.  2)  CT  3  =  C TV 
I F ( EPSLN  3  .EO.  888.  .ANT.  53  .FO.  898  .)  CT 3  =  C*PLX(0.0, 
GOTO  821 

RAY  IS  CLEAR  OF  THE  MASK  F2 
CONTINUE 

PSI1  =  DATANC  (Z1  -Z2)/  DC  ) 

CALL  FRESNL (EPSLN 1  ,  R L AMD  A , S 1 . PS  I , CTH , CT V ) 

CTH  =  CTH  *  REFL1 
CT V  =  CT V  *  REEL  1 
IF ( IPOL  .  EO  .  1  )  CT 1  =  CTH 
IF(IPOL  .  EO  .  2  )  CT 1  =  CTV 

IF( EPSLN 1  .EO.  888.  .AND.  SI  .EO.  888  .)  CT  1  =  CMPLXfO.O 
CONTINUE 

CALCULATE  CE82  USING  V82  CORRESP.  TO  F2 
CALL  DCS ( VC82 «  VS82,  V82) 

A  8  2  =  1 . /SORT  2  *  DSQPT ( ( VC  8  2  +0.5)**2  +CVS82  +0.5)**2) 

IF  ( VC8  2  +  0.5  .GE.  0.)  THEN 

R  =  DATANC  (VS82  +0.5)/(VC82  +0.5))  -PI/4. 

ELSE 

B  =  PI  +OATAN ( ( VS82  +0.5)/(VC82  +0.5))  -PI/4. 

END  IF 

A 1  =  A  8  2  *  A81  *  DSIN(B) 

B 1  =  A  8  2  *  A  8 1  *  OCOS(B) 

CE82  =  CMPLX(El.Al)  *  CT1  ♦  CT2  ♦  CT 3  ♦  GRdB 

CALCULATE  COMPLEX  F  (CF) 

PY  FIRST  CALCULATING  THE  COMPLEX  PARTS  CE11,  CF12, 

CF21,  CE31,  CE41,  CE52,  CE62,  CE71,  CE72,  CE81,  CE82 
CE32,  CE42.  CE22,  CE51,  CE61 

NOW  ADD  TO  GET  CF 

CF  =  +CE 1  +CE22  +CE32  +CF42  +CF52  +  CF.62  +  CF72  +CF92 
GOTO  216 

CF  =  CE11  +CE21  +CE31  +CF41 
WRITE  OUT  ANSWERS 

FO  =  CABS ( CF  ) 

FLOG  =  20.  *  A  LOG  1 0 ( F  0 ) 

WRITEI2,*) 

WRTTEC2,*) 

WRIT  E(2, 1 1 ) 

WRITE(2,12) 

WRITE (2, 14)  FLOG,  Z2 


WRITEI2, 19) 
WRITE(2 , 20) 
WR ITE ( 2 , 37  ) 


WRITE(2,18) 

VU  , 

V  2  1  , 

V  3  1  , 

V  4  1 

WRITE(2,38) 
WRITE(2, 18) 

V  1  2  , 

V22  , 

V  3  2  , 

V42 

WRITE ( 2 , 40  ) 
WRITEI2, 18) 

V  5  1  , 

V52, 

V  6  1  , 

V  6  2 

WRITE! 2 , 46 ) 
WRITE ( 2 , 1 8 ) 

V  7  1  , 

V7  2  , 

V  8  1  , 

VP2 

0.0) 


,0.0) 
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go  TO  100 

r  COMF  HERE  WHEN  FINISHEC. 

1000  CONTINUE 

WRITF(2,600)  Z2S,  Z2,  CZ2 
WR ITE ( 2 , 20 ) 

300  CONTINUE 
16  CLOSE ( 1 ) 

CLOSE ( 2 ) 

STOP 

1  FORMAT ( *  INPUT  Z1  ,  Dl,  D2,  03,  FI,  F2,  Z2S ,  DZ2'/) 

2  FORMAT  C '  INPUT  POLARIZATION..  H  =  HOR I ZON T A L , Vs V FR 1 1 C AL V  ) 

3  FORMAT ( A 1 ) 

4  FORMAT (  1  1  A  1  ) 


10  FORMAT ( 

ANTENNA  HEIGHT  ( M ) 

• 

• 

'  F6 . 2/ 

1 

OIST.  FROM  TRANSMITTER  TO  EOGE  FI  (M): 

'.FtO.l/ 

2 

DIST.  FROM  FI  TO  EDGE  F2 

OO : 

'  .F10.1/ 

3 

DIST.  FROM  F2  1C  TARGET 

RECEIVER  (M): 

',F10.1/ 

4 

HEIGHT  OF  KNIFE  EDGE  FI 

(M) : 

'.F10.1/ 

5 

HEIGHT  OF  KNIFE  EDGE  F2 

(M) : 

'  .F10.1/ 

6 

WAVELENGTH  (M)S 

'  ,F10.5/ 

7 

POLARIZATION: 

'  ,  A  1  / 

fl 

EPSILON  1  (EPSLN1) 

OVER 

Dl : 

'  ,  F8 . 3/ 

<3 

EPSILON  2  ( EPSLN 2  ) 

OVER 

D2 : 

'  ,F8.3/ 

1 

EPSILON  3  ( EPSLN 3  ) 

OVER 

D3 : 

'  , F8 . 3/ 

2 

SIGMA  1  (SI) 

OVER 

Dl : 

'  ,F10.3/ 

3 

SIGMA  2  (S 2) 

OVER 

D2: 

'  ,F10. 3/ 

4 

SIGMA  3  (S3) 

OVER 

D3: 

'  ,F10.3/ 

S 

* 

REFLECTION  (REEL  1  ) 

OVER 

Dl : 

',F10.3/ 

6 

• 

REFLECTION  (REFL2) 

OVER 

D2: 

-  .  F  1  0 . 3  / 

7 

• 

REFLECTION  (REFL3) 

OVER 

D3: 

'.F10.3//) 

600  FORMAT  C ‘ 

TARGET  HT  (M)  FRCM 

'.F5.C 

TO  ' , F8 . 4  ,  ' 

BY  '  , F5 , 1 / 

1  1 

FORMAT ( 55H 

POWER  GAIN  IN 

TARGET  HT 

POWER  i 

GAIN  OF  SCATTFRFP) 

12 

FORMAT (55H 

db,  20  x  LOG ( F )  (METERS) 

FIELD  9 

GROUND  ZERO  (db)) 

14 

FORMAT ( 3F1 5 

.4//) 

1  5 

F0RMAT(55HMAGNITUDES  of 

ELECTRIC  FIFLD 

COMPONENTS  ASSOCIATED  WITH) 

61 

FORMAT ( 41 H 

EACH  RAY  PATH 

,  EXPRESSED  IN 

volts/meter/) 

17 

FORMAT(55H 

E 1 1  MAG 

E21MAG 

E  3 1 maG 

E4  1  MAG  ) 

18 

FORMAT (4F1 5 

.4//) 

19 

FORMAT ( 55H 

INTEGRATION 

FARA“.  OF  FRESNFl  INTFGR ALS  DERIVED  FRCM) 

20 

FORMAT (25H 

Zl ,  Z2,  Dl  , 

C 2 ,  AND  D3/  ) 

22 

FORMAT (55H 

GRAZING  ANGLE  GRAZING  ANGLE 

GAIN  REDUCTION) 

23 

FORMAT (47H 

AT  SOURCE  (RAD)  AT  TARGET  (RAD) 

FACTOR/) 

25 

FORMAT ( 39H 

FRESNEL  REFLECTION  COEFFICIENT  OVER 

: ) 

26 

FORMAT (46H 

Dl 

D2 

D  3  / ) 

28 

FORMAT ( 35H 

PHASE  LAG  ANGLE  IN  DEGREES  OVER:) 

29 

FORMAT (2F1 5 

.4//) 

30 

FORMAT ( 5 1 H 

E12MAG 

E22MAG 

F.  32 MAG 

E42MAG) 

32 

F0RMAT(51H 

E  5  1  M  A  G 

E52MAG 

E61 MAG 

E62MAG) 

33 

FORMAT ( 27H 

Dl 

D2  /  ) 

34 

FORMAT ( 54H 

OFF 

-AXIS  BEAM  ANGLE) 

35 

FORMAT ( 49H 

TTHETA1  (RAC)  TTHETA2  (PAP)  OF 

ANTENNA  PHI) 

37 

FORM AT ( 49H 

VI  1 

V21 

V  3  1 

V  4 1  ) 

38 

FORMAT ( 49H 

V  1 2 

V22 

V  3  2 

V42 ) 

40 

FORMAT (49H 

V51 

V52 

V61 

V6?) 

41 

FORM AT (3F 15. 4/ ///////) 

42 

FORMAT ( 39H 

GRAZE  ANGLE 

GRAZE  ANGLE 

GRAZE  ANGLE) 

43 

FORMAT ( 55H 

AT  SOURCE 

CVER  D2 

OVER  D 3 

GAIN  REDUCTION) 
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44 

FORMAT ( 48H 

-RADI ANS- 

-KADI  ANS- 

-RAD  I ANS- 

F  ACTOR/ ) 

45 

FORMAT ( 5 1 H 

E71MAG 

E72MAG 

F81MAG 

F  R  2  m  A  G  ) 

46 

FORMAT ( 49H 

V  71 

V  72 

Vfi  t 

V  R  2  ) 

47 

FORMAT ( 8H 

F1MAG) 

48 

FORMAT (FI  5  . 

4//  ) 

FND 

************ 

************ 

************* 

************* 

************* 

C  SUBROUTINE  FOR  COMPLEX  REFLECTION  COEFFICIENTS 

SUBROUTINE  FRESNL (  El  ,  WAVE,  CONDUC,  ANG,  CTH ,  CTV) 
COMPLEX  CAK,  CTV,  CTV1,  CTV2,  CTH.  CTH1,  CTH2 


0  FI  .  THE  DIELECTRIC  CONSTANT  (FROM  0  TO  100) 

C  LAMDA  .  THE  WAVELENGTH  IN  METERS 

C  ,  CONDUC  .  THE  CONDUCTIVITY  IN  MHOS/METER 

C  ANG  .  THE  ANGLE  IN  RADIANS 


C  CALCULATE  THE  COMPLEX  CONSTANT 

AKI  =  -60.  *  WAVE  *  CONDUC 
CAK  =  CMPLX (£ 1 ,  AKI) 

C  CALCULATE  THE  VERTICAL  POLARIZATION  REFLECTION  COEFFICIENT 

CTV1  a  CAK  *  SIN ( ANG  )  -CSQRT(CAK  -COS ( A NG ) * *2  ) 

CTV2  a  CAL  *  SIN(ANG)  +CSORTCCAK  -COS ( ANG ) **2  ) 

CTV  =  CTV1  /  CTV2 

C  CALCULATE  THE  HORIZONTAL  POLARIZATION  REFLECTION  COEFFICIENT 

CTH 1  =  SIN(ANG)  -CSORTCCAK  -COS ( ANG ) *  *  2 ) 

CTH2  a  SIN (ANG)  +CSQRT(CAK  -COS ( A NG ) * * 2 ) 

CTH  a  CTH 1  /  CTH2 

1000  RETURN 
END 

C  SUBROUTINE  TO  EVALUATE  THE  FRESNEL  INTEGRALS 

SUBRCUTINE  DCS(C,  S,  X) 

IMPLICIT  REAL'S  (A-H,  0-Z) 

DIMFNSION  A  A  (  1 2  )  ,  RR(12),  CC(12),  DD(12) 

V  =  X 

PIE2  =  1 .57079632680  +  00 
X  =  PIF?  *  X  *  X 
7  a  DABS(X) 

I  F ( Z  . N  E .  0.0  +  00)  GOTO  1000 

C  =  O.D+OO 
S  =  O.D+OO 
X  a  U 
RETURN 

1000  CONTINUE 

I F ( Z  -  4.0D+00)  3,  3,  4 
3  C  a  DCOS(Z) 

S  =  DSIN(Z) 

Z  a  Z  /  4.00+00 
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PZ  =  DSORT(Z) 

ASUM  =  A  A ( 1 ) 

PSUM  =  BP ( 1 ) 

00  40  J  =  2.  12 

ASUM  =  ASUM  +  A  A  ( J )  *  Z**(J  -1) 

BSUM  =  BSUM  +  BBC J )  *  Z»*(J  -ID 
40  CONTINUF 

FC  =  DZ  ♦  (  S4BSUM  +C* ASUM  D 
FS  =  OZ  ♦  (»C*BSUM  +S*ASUM) 

C  =  FC 
S  =  FS 
GOTO  5 

4  CONTINUE 

P  =  DCOSCZ) 

S  =  DSIN(Z) 

Z  s  4.0D+00  /  Z 

csum  =  cccn 
psum  =  DD(i) 

DO  30  J  =  2,  12 

CSUM  =  CSUM  +CCCJ)  *  Z**(J  -1) 

30  CONTINUE 

PSUM  =  DSUM  +  DD  ( J )  *  Z**CJ  -1) 

Z  =  DSORT(Z) 

C  =  0.5D0  +Z  *  (0  *  CSU*  +  S  *  DSUM) 

S  =  0.5D0  +  Z  *  (S  *  CSUM  -C  ♦  DSUM) 

5  CONTINUE 
X  =  U 

I F ( U  .GT.  0  .  DO )  GOTO  6 
C  =  -C 

s  =  -s 

6  CONTINUE 
RETURN 

DATA  A  A  (  1 )  .  A A ( 2  )  /  0 . 1 5 9 5 76 9 1 4D+0 1  ,-0.1702D-05  / 

PATA  A  A  ( 3 )  ,  AAC4D  /-0 . 6 808568854D+O 1  ,-0.5763610-0  3  / 

PATA  A A  C 5 )  ,  A A  C 6  )  /  0 . 69 2 069 1 902D+0 1  , -0 . 1 6898657D-01  / 

DATA  A  A  (  7  )  ,  A  A ( ft  )  /-0 . 30 5 0 4 8 566D+0 1  , -0 . 7 5 7 5 2 4  1  9D -0 1  / 

PATA  AAC9)  ,  A  A (  1 0  ) /  0 . 8 5 0 6 6 3 7 8 1 D +  00  ,-0.256390410-01  / 

DATA  A  A  ( 1 1 ) ,  AAC12D/-0.15023096D  +  00  ,  0.344047790-01  / 

PATA  BO(l)  ,  80  C  2 )  /-0. 330-07  ,  0.42553875240  +  01/ 

PATA  BB ( 3 )  ,  BB ( 4 )  /-0. 928  10-04  , -0. 77800204D  +  01  / 

PATA  BB  (5 )  ,  PBC6)  /-O. 95208950-0?  ,  0 . 507 5 1 6 1 29  ft 0  +  0 1 / 

PATA  BB ( 7 )  ,  BB  (  8  )  /-0 . 1  3 8 3 4 1  9 4 7 D  +  0 0  , - 0 . 1  3 6 3 7 2 9  1  2 4D  +  0 1 / 

PATA  BB (9 )  .  8BC10D/-0.403349276D+00  ,  0 . 70222 20 1 6D+00  / 

PATA  BB(ll),  BB(  12  )/-0.216  195929D  +  00  ,  0 . 1  954703  1  C-Ot  / 

DATA  CC(l)  .  CC ( 2 )  /  0.00+00  ,-0.249339570-01  / 

PATA  CC ( 3 )  ,  CCC4)  /  0.3936D-05  ,  0 . 57 70956D-02  / 

DATA  CCC5D  ,  CC ( 6  )  /  0.689892P-03  ,-0.94971  360-02  / 

PATA  CC (7 )  ,  CC(fl)  /  0.119488090-01  ,-0.67488730-02  / 

PATA  CC (9 )  ,  CCC10)/  0.246420-03  ,  0.21029670-02  / 

DATA  CC(ll),  CCC12D/-0.121793D-02  ,  0.233939D-03  / 

PATA  DDC1)  ,  DD ( 2  )  /  0 . 1 994 7 1  1  4P  +  00  ,  0.23D-07  / 

PATA  DOC  3 )  ,  DDC4)  /-0 . 9 35 1  3 4 1 D-02  ,  0.23006C-04  / 

PATA  DD ( 5 )  ,  OD ( 6 )  /  0.4851  46D-02  ,  0.1  90321  80-02  / 

DATA  DD ( 7 )  ,  DD(8  )  /-0 . 1 7 1 229 1 4D-0 1  ,  0.290640670-01  / 

DATA  DDC9)  ,  00 C 1 0 ) /-0 . 2 7928955D-0 1  ,  0.164973080-01  / 

DATA  DD ( 11),  DD(  12D/-0. 55985150-02  ,  0.fl383fl6D-03  / 


Appendix  C. 


Computer  Program  DOUBLE_PLOT 
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r  PROGRAM  DOUBLE- PLOT 

CHARACTERS  1  DUMMY 

C  PROGRAM  TO  SET  UP  HT  VS.  dB  VALUES 

OPEN(l,  NAME  =  'DOUBLEDGE.OUT' 
OPEN (  2 ,  NAME  =  ' PLOT . DAT  ' 

D  =  3001 

DO  1  I  =  1,  0 

DO  3  J  *  1|  23 

READ (1,99)  DUMMY 

3  CONTINUE 

4  READ ( 1 , * ,  END  =  2)  FLOG 

WRITE(2,*)  FLOG 

1  CONTINUE 

2  CLOSE ( 1 ) 

CL  OSE ( 2  ) 

STOP 

Q9  FORMAT ( A1  ) 

FND 


FOP  EZGRAPH 


TYPE  =  'OLD') 
TYPE  s  'NEW '  ) 


,  Z2 
,  Z2 
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Appendix  D.  —  Plots  of  Computer  Studies  of  Doubly  Diffracted  Horizontally 
Polarized  Microwaves 


This  appendix  presents  plots  of  horizontally  polarized  microwaves  with 
the  hill  closer  to  the  source  set  at  a  15-m  height  and  the  hill  farther  away 
from  the  source  varying  in  height  from  5  to  45  m. 

These  figures  show  the  effect  of  terrain  on  a  source  with  nonuniform  an¬ 
tenna  gain.  Gain  of  power  density  relative  to  free  space  is  shown  in 
decibel  scale  versus  target  height 
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STUDY  OF  POWER  LOBE  STRUCTURE 


Figure  D-2.  Case  of 
hill  IS  m  high  near 
source  and  hill  10  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigmo=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflecfion 
g  hill  hts:  0=15  m,R=10  m.  freq=10  GHz,  Dl=2500.  D2=2500.  D3=2500,ont  ht=3  m,  gain 


APPENDIX  D 


STUDY  OF  POWER  LOBE  STRUCTURE 


moist  ground:  sigma=3.,  dielec.  const.=10.,  "H"  polarization;  diffuse  reflection 


Figure  D-4.  Case  of 
hill  15  m  high  near 
source  and  hill  20  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigma=3.,  dielec.  consf.=lO.,  "H"  polarization;  diffuse  reflection 
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STUDY  OF  POWER  LOBE  STRUCTURE 


Figure  D-5.  Case  of 
hill  15  m  high  near 
source  and  hill  25  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


moist  ground:  sigma=3.,  dielec.  const.=10.,  "H"  polarization;  diffuse  reflection 
9  hill  hts;  1=15  m.  R=25m,  freq=10  GHz,  01=2500,  02=2500,  D3=2500,ant  ht=3  m,  gain 


-60.0  -50.0 


-J0.0  -20.0 

db 


Figure  D-6.  Case  of 
hill  15  m  high  near 
source  and  hill  30  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


moist  ground:  sigmo=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
9  hill  hts:  L=15m,  R=30m.  freq=10  GHz.  D1=2500.  D2=2500,  D3=2500.ont  ht=3  m.  goin 
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STUDY  OF  POWER  LOBE  STRUCTURE 


Figure  D-7.  Case  of 
hill  15  m  high  near 
source  and  hill  35  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


moist  ground:  sigma=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
9  hill  hts:  1=15  m.  R=35m,  freq=10  GHz.  01=2500,  02=2500,  D3=2500.ant  ht=3  m,  gain 


Figure  D-d.  Case  of 
hill  15  m  high  near 
source  and  hill  40  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigmo=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
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STUDY  OF  POWER  L08E  STRUCTURE 


Figure  D-9.  Case  of 
bill  15  m  high  sear 
source  and  hill  45  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is3m. 


moist  ground:  sigma=3„  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
hill  hts:  1^15  m,  R=45m,  freq=10  GHz,  Dl=2500,  D2=2500,  D3=2500.ont  nt=3  m,  gain 


-70.0 


-60.0 


-50.0 
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-50.0 
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-10.0 
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dB 


Figure  D-10.  Case  of 
hill  30  m  high  near 
source  and  hill  5  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigma=3.,  dielec.  const.=10.,  "H"  polarization;  diffuse  reflection 
9  hill  hts:  U=30  m,  R=  5m,  freq=10  GHz,  D1=2500.  D2=2500,  D3=2500,ant  ht=3  m,  gain 
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STUDY  OF  POWER  LOBE  STRUCTURE 


Figure  D-12.  Case  of 
hill  30  m  high  near 
source  and  hill  15  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigmo=3.,  dielec,  const. =10.,  "H"  polorization;  diffuse  reflection 
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STUDY  or  POWER  LOBE  STRUCTURE 


Figure  D-13.  Owe  of 
hill  30  m  high  near 
source  and  hill  20  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


moist  ground:  sigma=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
hill  hts:  L=30  m,  R=20m,  freq=10  GHz,  01=2500,  D2=2500,  D3=2500,an1  ht=3  m,  gain 


-60.0  -50.0 


Figure  D- 14.  Case  of 
hill  30  m  high  near 
source  and  hill  25  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


moist  ground:  sigma=3.,  dielec.  const.=10.,  "H"  polarization;  diffuse  reflection 
°  hill  hts:  1=30  m,  R=25m,  freq=10  GHz,  01=2500,  D2=2500.  D3=2500,ant  ht=3  m.  gain 
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Figure  D-15.  Case  of 
hill  30  m  high  near 
source  and  hill  30  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


STUDY  OF  POWER  LOBE  STRUCTURE 

moist  ground:  sigmo=3.,  dielec,  const. =10.,  "H"  polarization;  diffuse  reflection 
9  hill  hts:  L=30  m,  R=30m,  freq=10  GHz,  01=2500,  02=2500,  D3=2500,ant  ht=3  m,  gain 


Figure  D-16.  Case  of 
hill  30  m  high  near 
source  and  hill  35  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


—37.5  -25.0 

dB 


moist  ground:  sigma=3.,  dielec.  const.=10.,  ”H"  polarization;  diffuse  reflection 
9  hill  hts:  L=30  m,  R=35m,  freq=10  GHz.  01=2500,  D2-2500,  D3=2500,ont  ht=3  m.  gain 
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STUDY  OF  POWER  LOBE  STRUCTURE 


Figure  D- 17.  Case  of 
hill  30  m  high  near 
source  and  hill  40  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 
is  3  m. 


moist  ground:  sigma=3.,  dietec.  const. =10.,  "H"  polarization;  ditfuse  reflection 

o  hill  hts:  1=30  m,  R=40m,  freq=10  GHz,  01=2500,  02=2500,  03=2500,ant  ht=3  m,  gain 
<?-, - . - - — _ - - - - — - 1 


-S7.S  -25.0 

dB 


moist  ground:  sigma=3.,  diefec.  const.=t0.,  "H"  polarization;  diffuse  reflection 
hill  hts:  1=30  m,  R=45m,  freq=10  GHz,  01=2500,  D2=2500,  D3=2500,ant  ht=3  m.  gain 


Figure  D- 18.  Case  of 
hill  30  m  high  near 
source  and  hill  45  m 
high  near  target,  with 
horizontal  polariza¬ 
tion.  Antenna  height 


Appendix  E.  —  Power  Density  Comparison  Plots  at  Target: 

(1)  Double  Diffraction  Versus  Single  Diffraction  (fig.  E-l  to  E-3) 

and 

(2)  Transmitter  Height  of  3  m  Versus  6m  (fig.  E-4  to  E-7) 
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Figure  E-l.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over 
two  edges  to  power 
density  loss  resulting 
from  diffraction  over 
one  edge.  Positive  or 
negative  gain  in  deci¬ 
bel  units  is  measured 
relative  to  free  space 
propagation  at  0  dB. 
In  this  case,  single 
edge  is  set  closer  to 
microwave  emitter 
than  to  target. 


Figure  E-2.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over 
two  edges  to  power 
density  loss  resulting 
from  diffraction  over 
one  edge.  Positive  or 
negative  gain  in  deci¬ 
bel  units  is  measured 
relative  to  free  space 
propagation  at  0  dB. 
In  this  case,  single 
edge  is  set  midway 
between  emitter  and 
target 


POWER  LOSS  OF  DOUBLE  VS.  SINGLE  DIFFRACTION 
conduc.=3  s,  permiH.=10,  "H”  poloriz.,  reflect. =.3,  hills=15m,  ont.  ht.=3m 


db 


db 
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Figure  E-3.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over 
two  edges  to  power 
density  loss  resulting 
from  diffraction  over 
one  edge.  Positive  or 
negative  gain  in  deci¬ 
bel  units  is  measured 
relative  to  free  space 
propagation  at  0  dB. 
In  this  case,  single 
edge  is  set  closer  to 
target  than  to 
microwave  emitter. 


Figure  E-4.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over 
two  I5-m-high  edges 
for  microwave  emitters 
set  at  heights  of  3  and 
6  m,  respectively. 

Solid  curve  shows  an¬ 
tenna  height  at  3  m, 
and  dotted  curve,  at 
6  m. 


POWER  LOSS  OF  DOUBLE  VS.  SINGLE  DIFFRACTION 
conduc.=3  s.  permitt.=10,  "H"  poloriz.,  reflect. =.3,  hills=15m,  ant.  ht.=3m 


-60.0  -50.0  -40.0  -30.0  -20.0  -10.0  0.0  >0.0 

db 

POWER  LOSS  OF  3  m.  vs.  6  m.  ANTENNA  HEIGHT 


conduc.=3  s,  permitt.=10,  "H"  poloriz.,  reflect.=.3,  hills=15  m 


-70.0  -60.0  -30.0  -40.0  -30.0  -20.0  -10.0  0.0  K.O 

db 
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Figure  E-5.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over  a 
single  knife  edge  set  at 
a  15-m  height  to 
microwave  emitters  set 
at  3-  and  6-m  heights, 
respectively.  Solid 
curve  shows  antenna 
height  at  3  m,  and 
dotted  curve,  at  6  m. 
Edge  is  here  set  closer 
to  emitters  than  it  is  to 
target 


Figure  E-d.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over  a 
single  knife  edge  set  at 
a  15-m  height  to 
microwave  emitters  set 
at  3-  and  6-m  heights, 
respectively.  Solid 
curve  shows  antenna 
height  at  3  m,  and 
dotted  curve,  at  6  m. 
Edge  is  here  set  mid¬ 
way  between  emitters 
and  target. 


POWER  LOSS  OF  3  m.  vs.  6  m.  ANTENNA  HEIGHT 
co nduc.=3  s,  permitt.=10,  "H“  poloriz.,  reflect. =.3,  hi!ls=15m 


conduc.=3  s,  permiti.=10,  "H"  poloriz.,  reflect.=.3,  hills=15m 
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POWER  LOSS  OF  3  m.  vs.  6  m.  ANTENNA  HEIGHT 


Figure  E-7.  Compari¬ 
son  (over  electrically 
interactive  terrain)  of 
power  density  loss 
from  diffraction  over 
a  single  knife  edge  set 
at  a  15-m  height  to 
microwave  emitters  set 
at  3-  and  6-m  heights, 
respectively.  Solid 
curve  shows  antenna 
height  at  3  m,  and 
dotted  curve,  at  6  m. 
Edge  is  here  set  closer 
to  target  than  to 
microwave  emitters. 


conduc.=3  s,  permift.=lO,  "H"  poloriz.,  reflect.=.3,  hills=15m 
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Appendix  F.  —  Plot  Comparing  a  Theoretically  Derived  Curve  for  Two  Dif¬ 
fractions  over  Terrain  with  a  Curve  Prepared  from  Field  Measurement  Data 
for  Two  Diffractions  over  Terrain 
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Figure  F-l.  Compari¬ 
son  of  theoretically  de¬ 
rived  curve  for  two 
consecutive  diffrac¬ 
tions  with  curve  de¬ 
rived  from  field  meas¬ 
urement  data  for  sig¬ 
nal  transmission  over 
two  peaks.  Path  is 
Forest  Hill  azimuth 
relative  to  transmitter. 
Transmitter  is 
992.95  m  above  sea 
level;  ground  zero  is 
276  m  above  sea  level. 


FOREST  HILL  AZIMUTHAL  PATH  ALONG  TERRAIN 


Ft  112m,  F2:  72m,  anl.  &  116.95m  above  ground  zero,  lambda  =  2.7125m 

av.  of  VHF  meas.  1,  2  0  5m  inlervals;  H  polariz.,  D  1=3.4  tkm,D2=3.56km,  D3=3.62km 
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Appendix  G.  —  Plot  Comparing  a  Theoretically  Derived  Curve  for  a  Single 
Diffraction  over  Terrain  (using  the  KNIFEDGE  Program)  with  a  Curve 
Prepared  from  Field  Measurement  Data  for  Two  Diffractions  over  Terrain 
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APPENDIX  G 


Figure  G-l.  Compari¬ 
son  of  theoretically  de¬ 
rived  curve  for  a  single 
diffraction  with  curve 
derived  from  field 
measurement  data  for 
signal  transmission 
over  two  peaks.  Path 
is  Forest  Hill  azimuth 
relative  to  transmitter. 
Transmitter  is  392.95 
m  above  sea  level; 
ground  zero  is  276  m 
above  sea  level. 


STUDY  OF  POWER  LOBE  STRUCTURE  (Forest  Hill  azimuth) 


Comparsion  of  M.l.T.  double  diffraction  doto  with  single  edge  theory 
sigmo=.005  mho/m,  dielec.  const.=15.,  "H"  polarization,  lambda  =  2.7125 
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